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ABSTBACT 


A  direct  design  method  for  solving  the  problem  of 
robustness  to  cross-coupling  perturbations  in  multivariable 
control  systems  is  presented.  The  method  uses  numerical 
optimization  procedures  to  manipulate  the  system  feedback 
gains  as  direct  design  variables.  The  manipulation  is  accom¬ 
plished  in  a  manner  that  produces  desired  performance  by 
pole  placement  and  robustness  by  modification  of  the  minimum 
singular  values  of  the  system  return  difference  matrix. 

Channels  affected  by  cross-coupling  perturbation  may  be 
recognized  by  the  character  of  their  transfer  function  plot. 
The  mechanism  used  by  the  pole  placement  and  robustness 
routine  in  obtaining  a  robust  design  is  evident  from  the 
gain  changes  associated  with  the  transfer  function  diagram 
and  the  zero  shifts  shown  cn  pole-zero  plots.  The  pole 
placement  and  robustness  routine  uses  gain  equalization  and 
zero  assignment  to  modify  the  characteristics  of  the  system 
in  the  areas  of  low  singular  values,  producing  a  robust 
design. 

A  modification  of  the  pole  placement  and  robustness 
routine  that  may  be  applied  to  the  design  of  robust 
observers  is  also  presented.  Using  feedback  and  filter  gains 
as  direct  design  variables  a  practical  design  procedure  for 
robustness  recovery  in  observer  based  systems  is  obtained. 
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I.  18  TRODUCTIQN 


With  the  rising  interest  in  multivariable  control  theory 
brought  on  by  increasingly  complex  systems  the  need  has 
arisen  to  develop  design  methods  that  will  allow  the 
designer  tc  specify  system  performance  while  at  the  same 
time  ensuring  relatively  high  stability  margins  or  robust¬ 
ness.  In  the  single-input  single-output  (SISO)  case  the 
designer  has  had  the  tools  to  do  these  tradeoffs  in  the  form 
of  Nyguist,  Bode  and  root  locus  plots.  In  the  multi-input 
multi-output  (MIMO)  case  the  classical  methods  are  not 
totally  appropriate. 

With  the  increased  interest  in  MIMO  systems  numerous 
methods  of  design  have  been  employed  to  obtain  suitable 
system  performance  and  robustness  with  varying  degrees  of 
success.  One  primary  method  of  design  is  to  keep  the  plant 
as  decoupled  as  possible  throughout  the  design  so  that  each 
individual  element  may  be  controlled  •  independently  and 
designed  essentially  as  a  single  loop  system.  Rosenbrock 
[Hef .  1  ]  has  developed  a  procedure  where  the  multiloop 
system  is  modified  into  a  system  that  has  diagonal 
elements  that  are  much  larger  than  any  off-diagonal 
elements.  This  diagonally  dominant  system  is  then  in  a  form 
where  conventional  Nyguist  type  techniques  can  be  employed 
in  the  analysis.  A  third  common  HlflO  design  method  is  that 
of  the  linear  Quadratic  (L  Q)  method.  This  method  uses  a 
quadratic  cost  functional  and  optimization  principles  to 
allow  the  designer  tc  design  for  various  performance  levels 
by  adjusting  the  matrix  weighting  terms  used  in  the  cost 
function.  The  major  difficulty  with  all  of  the  above 
methods  is  that  they  are  not  necessarily  robust .  This  is 
especially  true  for  cross-coupling  terms  between  loops. 


The  primary  achievement  of  this  thesis  has  been  the 
incorporation  of  the  time  domain  pole  placement  design 
procedure  with  a  method  of  using  return  difference  matrix 
singular  values  to  improve  the  robustness  of  the  design.  Ihe 
technigue,  which  utilizes  a  modern  optimization  routine,  can 
significantly  assist  the  designer  in  obtaining  robustness  in 
the  face  of  cross-coupling  perturbations.  It  has  also  been 
shown  that  the  cross-coupling  perturbation  problem  can  be 
detected  by  using  classical  open-loop  •  ?de  diagrams  as  well 
as  modern  control  analysis.  The  pole  placement  and  robust¬ 
ness  design  routine  developed  for  this  thesis  has  been  used 
cn  several  problems  discussed  in  recent  literature.  In  these 
studies  the  pole  placement  and  robustness  design  code  has 
proven  capable  of  meeting  the  desired  goals  cf  pole 
placement  and  robustness  and  also  brought  to  light  seme 
interesting  aspects  of  the  cross-coupling  perturbation 
problem.  A  slightly  modified  pole  placement  and  robustness 
routine  has  proven  effective  in  the  design  of  robust 
observers. 

The  remainder  of  the  thesis  will  present  background 
material  on  SISO  systems  in  Chapter  Two  and  on  HIMO  systems 
in  Chapter  Three.  Optimization  will  be  discussed  in  Chapter 
Four  along  with  a  discussion  of  the  Automated  Design 
Synthesis  (ADS)  program  used  as  the  optimizer  routine  for 
the  pc]e  placement  and  robustness  design  procedure  developed 
in  this  thesis.  The  thesis  methodology  will  be  discussed  and 
outlined  in  Chapter  Five.  This  will  be  followed  by  chapters 
discussing  results  ficm  selected  problems.  Conclusions  will 
be  presented  in  the  final  chapter. 


II.  SINGLE-INPUT  SINGLE-OUTPUT  SYSTEMS 


The  purpose  of  this  thesis  research  has  been  to  develop 
a  method  of  obtaining  a  robust  multivariable  control  system 
design.  &  brief  review  of  the  concept  of  robustness  and 
stability  in  the  framework  of  a  conventional  SISO  system 
will  be  done  before  pursuing  the  concepts  in  a  more 
complicated  fashion  in  the  following  chapters.  A  simple 
interpretation  of  robustness  is  the  ability  of  the  system  to 
tolerate  design  perturbations.  These  perturbations  could  be 
in  the  form  of  actuator  failures,  plant  parameter  uncer¬ 
tainty,  unmodeled  dynamics  or  nonlinear  terms,  or  any  one  of 
many  other  perturbations  to  the  nominal  design  of  the 
system. 

The  primary  reason  for  feedback  systems- is  the  control 
cf  uncertainty  within  the  system.  By  appropriate  use  of 
feedback,  properties  that  would  lead  to  an  unstable  system 
may  be  controlled.  When  stability  and  robustness  aspects 
are  considered  for  a  SISO  system,  frequency  domain  design 
concepts,  using  either  Nyguist  or  Bode  plots,  are  normally 
used.  Robustness  in  SISO  systems  is  formulated  naturally  by 
the  concept  of  gain  and  phase  margins,  both  of  which  are 
readily  available  on  the  Nyguist  or  Bode  diagram. 

In  figure  2.1a  nominal  feedback  system  can  be  seen  with 
a  perturbation  element  (s)  placed  in  series  with  the 

nominal  system.  When  of  =1  the  system  is  nominal  and  stable. 
To  determine  the  positive  phase  margin  of  the  system  the 
value  of  (s)  =  c({  jw)  =  eJ<^  will  be  changed  oy  varying 
until  the  system  just  becomes  unstable.  This  value  of  if  will 
then  be  the  system  phase  margin.  The  negative  phase  margin 
can  be  computed  in  the  same  way.  To  find  the  gain  margin  the 
magnitude  value  of  is  increased  until  the  system  just 
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Figure  2.1  Basic  Gain  and  Phase  Kargin  Determination  Hodel. 

becomes  unstable.  This  is  the  upward  gain  margin.  A  similar 
lower  margin  is  also  defined. 

These  gain  and  phase  margins  may  not  be  adequate  meas¬ 
ures  cf  robustness  [Ref.  2]  because  they  do  not  account  for 
simultaneous  variation  in  both  gain  and  phase.  Therefore, 
while  large  individual  gain  or  phase  changes  may  not  desta¬ 
bilize  the  system,  small  simultaneous  changes  in  gain  and 
phase  may  destabilize  the  system.  This  is  not  a  major  diffi¬ 
culty  in  classical  SISO  techniques  because  the  effect  can  be 
easily  detected. 

Gain  and  phase  margin  can  be  defined  in  terms  of  the 
cpen-lcop  frequency  domain  plots  in  either  the  Bode  or 
Nyquist  format.  Figure  2.  2  depicts  a  classical  Bode  plot 
showing  gain  and  phase  margin  determination  from  the  plot. 
The  Nyquist  plot  may  also  be  used  to  obtain  this  informa¬ 
tion.  Nyquist  criterion  states  that  if  the  cpen-loop 
transfer  function  G  (s)  H  (s)  has  n  poles  in  the  right  half 
plane  and  the  limit  cf  G  (s)  H  (s)  =constant  as  s-*«o  then  for  a 
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Figure  2.2  Classical  Bode  Plot. 

stable  system  the  locus  of  G(s)H(s)  will  encircle  the  -1+jO 
point  n  times  in  the  counterclockwise  direction  as  s  varies 
along  the  Nyguist  contour.  If  there  are  no  pales  in  the 
right  half  s  plane  then  the  locus  will  not  encircle  the  -1  ♦ 
jO  point.  Ihe  diagram  in  figure  2.3  illustrates  a  nominally 
stable  system.  The  gain  and  phase  margin  may  be  determined 
directly  from  the  diagram. 

Any  change  in  the  loop  transfer  function,  provided  the 
order  of  G  (s)H (s)  does  not  change,  that  changes  the  number 
of  times  the  locus  of  G(s)H(s)  encircles  the  (-1,0)  point  in 
the  Nyguist  plot  causes  the  system  to  become  unstable.  This 
leads  to  the  conclusion  that  the  minimum  distance  of  the 


Figure  2.4  Additively  Perturbed  Systea. 

Since  the  system  is  stable  the  {-1,0}  point  is  encircled 
the  correct  number  of  times  by  the  nominal  plant.  If  the 
locus  of  g(j<A>)  in  the  diagram  is  warped  until  it  passes 
beyond  the  (-1,0)  point  then  clearly  the  number  of  encircle¬ 
ments  of  this  point  will  change  and  the  system  will  become 
unstable,  assuming  the  order  of  the  plant  is  not  changed  by 
the  perturbation.  To  keep  the  locus  of  points  from  moving 
beyond  the  (-1,0)  point  eguation  2.1  must  hold. 


18 


I  £g  (juj)  |<|  1  ♦  g  (joJ)  |  (2.1) 

This  condition  is  illustrated  in  figure  2.6.  The  right-hand 
side  of  equation  2.1  is  just  the  magnitude  of  the  return 
difference  transfer  function  of  the  nominal  system.  The 
multiplicative  case  is  depicted  in  figure  2.7  with  its  asso¬ 
ciated  Nyguist  plot  in  figure  2.8.  The  requirement  for 
stability  is  similar  to  the  additive  case  and  may  be  stated 
in  eguation  2.2 

1  <l9(jw)  I  <1  1  ♦  <g(jw)-»)  J  (2.2) 
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The  above  arguments  will  be  applied  again  in  Chapter  3  to 
develop  aultivariable  stability  and  robustness  properties. 

The  linear  quadratic  design  has  been  the  primary  method 
employed  in  modern  control  design  practice.  In  this  method 
an  optimal  state  feedback  control  lav  is  developed  to  find  a 
set  of  feedback  gains  that  optimizes  a  chosen  performance 
index.  The  performance  index  for  the  steady-state  case  is 
given  in  equation  2.3 
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U(S)  E(S) 


+ag(s) 


G(S) 


Y  (s) 


Figure  2.7  Multiplicative  System. 


PI  =  /(?'Q  ?  ♦  u*B  u)dt 


(2.3) 


where  the  x’x  tern  and  the  u'u  term  form  quadratics.  The  Q 
and  B  matrices  are  chosen  by  the  designer  to  provide  the 
best  compromise  between  the  minimum  error  of  the  system  and 
the  minimum  energy  needed  to  control  the  system.  The  LQ 
method  is  based  on  the  use  of  closed-loop  state  variable 
feedback  fcr  the  control  of  the  system.  In  the  MIMO  problem 
IQ  methods  have  been  used  extensively  because  of  their  guar¬ 
anteed  stability  margins  with  diagonal  weightiag  matrices. 
For  diagonal  weighting  matrices  the  LQ  method  yields  a 
guaranteed  phase  margin  of  60  degrees  and  -6  db  to  infinite 
gain  margin. 


Figure  2,8  Nyguist  Plot  for  Hultiplicative  System. 

For  SISO  linear  quadratic  optimal  regulators  these 
stability  margins  can  be  developed  from  inequality  2.4 
[Bef.  3] 

|  1  ♦  f »  (joJI  -  a)  -*ti>  1  (2.4) 
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Writing  the  Kalman  eguation  as  2.5 


M  ♦  f •  (sl-a)-ib|  z  =  1  ♦  (1/^)1  G^s)!*  (2.5) 

cne  has  further 

|  1  ♦  G (s ) H (s)  |  2  =  1  +  (1/^)  JGK{S)  |2  (2.6) 

and  for  all  s*j  u>  and  0<  <o  <  »  the  function  (V^)  G<  (s)  is 
greater  than  zero,  therefore  the  Kalman  inequality  is  shown 
to  be  2.7. 

|  1  +  G  { s )  H (s)  |  >  1  (2.7) 
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The  graphical  interpretation  of  this  result  is  that  the 
polar  plot  cf  G{s)H(s)  must  remain  outside  the  unit  circle 
centered  at  the  -1  +  jO  point  for  all  frequencies.  Figure 


Figure  2.  10  Polar  Plot  of  an  Optimal  System. 

2.10  shows  such  a  polar  plot.  Since  the  optimal  regulator 
with  the  loop  broken  at  the  input  to  the  plant  (denoted  by 
the  x  in  figure  2.9  )  does  not  penetrate  the  unit  disk  about 
-1  +  jO  this  means  that  the  single  input  regulator  will  have 


a  phase  margin  of  at  least  60  degrees  and  a  gain  margin 
tolerance  of  fifty  percent  gain  reduction  and  infinite 
upward  margin.  Further  discussion  of  this  property  may  be 
found  in  Anderson  and  Moore  £Hef.  4].  With  this  basic 
review  of  the  concepts  of  stability  and  robustness  in  the 
classical  SISO  system  complete,  the  next  chapter  will  extend 
some  of  these  basic  concepts  to  the  MI MO  system. 


III.  MOLTI VARIABLE  SYSTEMS 


linear  Quadratic  design  has  performed  relatively  well  in 
aircraft  control  concepts  because  of  the  ability  tc  formu¬ 
late  the  system  state  equation  and  to  quantify  acceptable 
performance  indices  for  the  system.  Industrial  applications 
for  IQ  theory  have  been  less  successful  and  have  led  British 
researchers  to  look  into  forms  of  decoupled  design  methods. 
Cne  of  the  primary  methods  used  in  multivariable  design  is 
to  make  the  system  totally  decoupled.  This  method  allows 
each  loop  tc  be  designed  as  a  separate  entity  by  classical 
means.  Cne  of  the  primary  difficulties  with  this  method  is 
the  problem  of  finding  a  compensator  which  will  totally 
decouple  the  system.  The  method  also  suffers  from  the 
effects  of  cross-coupled  perturbation  terms.  A  method  that 
does  not  totally  decouple  the  system  but  makes  the  design 
problem  simplier  by  designing  a  compensator  that  only  causes 
the  diagonal  terms  of  the  transfer  matrix  to  be  dominant 
over  all  cff-diagonal  terms  has  also  been  developed, 
[Ref.  5].  Classical  frequency  domain  techniques  are  then 
used  to  design  each  loop  of  the  system.  The  major  difficulty 
with  this  and  other  single  loop  design  techniques  is  their 
failure  to  account  for  cross-coupling  perturbation  terms 
that  may  interact  between  the  loops.  The  individual  loops 
may  be  highly  robust  in  these  designs  but  the  overall  system 
robustness  may  be  low  because  of  the  loop  interaction  not 
accounted  for  in  the  design.  This  is  the  precise  area  that 
the  singular  value  analysis  has  proven  so  beneficial  in 
design  .  Singular  value  concepts  may  be  applied  tc  conven¬ 
tionally  designed  systems  to  assess  their  robustness.  for 
instance,  if  a  system  is  designed  by  L2  methods  the  designer 
may  then  formulate  the  transfer  function  of  the  system  and 


assess  tie  singular  values  of  the  return  difference  matrix 
of  the  system.  If  the  minimum  singular  value  is  found  to  be 
low  at  some  critical  frequencies  the  designee  can  then 
modify  the  Q  and/or  E  matrices  chosen  in  the  LQ  performance 
index  and  recalculate  the  design.  In  this  iterative  fashion 
a  robust  design  would  be  developed. 

A  generalization  of  the  SISO  Nyquist  theory  discussed  in 
the  previous  chapter  has  been  made  for  the  MIHO  problem. 
This  generalization  leads  directly  to  the  singular  value 
concept.  The  generalization  is  expressed  in  the  form  of  the 
multivariable  Nyquist  theorem  which  requires  that  a  closed 
loop  stable  system  have  the  same  number  of  counterclockwise 
encirclements  of  the  origin  by  the  locus  of  the  det  (I+g  (jej) ) 
as  the  number  of  open  loop  poles  that  are  unstable.  This 
theoren  is  formally  stated  as; 

let  N[f  (s)  ]  denote  the  number  of  clockwise  encirclements  of 
(-1,0)  by  the  locus  of  f  (s )  as  s  traverses  the  contour  D 
of  figure  3.1  in  a  clockwise  sense.  The  closed-loop  system 
will  be  stable  if  and  only  if  for  all  R  sufficiently  large 
N[f(s)]  =  -P 

where 

f(s)  =  -1  ♦  det[  I+G  (s)  J  =  ^(s)/  C/’nts )  -1  and 

P  =  the  number  of  closed  right-half  plane  zeros  of 
(s)  . 

The  application  cf  the  Nyquist  theorem  comes  through  the 
fact  that  a  multivariable  system  will  not  be  robust  to 
modelling  errors  if  the  return  difference  matrix,  I  +  g,  is 
nearly  singular  for  some  frequency.  If  I  +  g  is  nearly 
singular  a  small  change  in  g  may  make  I  +  g  exactly 
singular.  This  causes  the  det  (I  +  G)  to  become  zero  and  the 
Nyquist  encirclement  count  to  change  indicating  an  unstable 
system.  It  is  possible  for  very  small  changes  in  I+G  to 
produce  large  changes  in  the  determinant  of  I+G.  The  matrix 
I+G 


X:  Zeros  <^L 


Figure  3.1  Myguist  D  Contour. 

[10/s+a  9.99/s+aj 
10/s*a  10/s+a  j 

has  determinant  0.  1/ (s+a)  2.  If  the  element  is  changed  by 
only  cne  percent  to  9. 9/s+a  the  determinant  becomes 
1.1/ (s+a)2  which  is  a  significant  change  in  the  determinant 
value.  Therefore,  it  is  evident  that  det  (I  ♦  G>  is  net  an 
accurate  measure  of  how  near  the  return  difference  is  to 
singularity.  Researchers,  in  the  field  of  controls 
£Ref.  6],  [Bef.  7],  [Ref.  8],  [Bef.  9]  have  used  singular 
value  analysis  to  determine  how  near  the  return  difference 
matrix  is  tc  singularity. 

Since  the  number  cf  encirclements  of  the  Nyguist  diagram 
changes  as  f  (s)  passes  through  the  -1  point  or  when  det(I+G) 
is  zero  it  is  important  to  find  how  near  the  return  differ¬ 
ence  matrix  I+G  is  to  being  singular.  This  nearness  to 
singularity  can  be  interpreted  as  closeness  of  the  matrix  G 


to  the  critical  point,  -1.  A  quantity  which  can  be  used  to 
express  the  nearness  to  singularity  of  the  matrix  is  the 
minimum  matrix  singular  value  denoted  by  (f  .  Given  a  matrix 
A  the  singular  value  may  be  expressed  by  equation  3. 1 

C  (A)  =  min  (  (|-A))"*‘  (3.1) 

where  \  is  the  eigenvalue  of  the  complex  conjugate 

transpose  of  A  times  A.  A  basic  MIMO  linear  system  is 


Figure  3.3  Additive  Perturbation. 

encircle  the  -1  point  the  system  will  remain  stable.  A 
sufficient  condition,  recalling  the  SISO  discussion  in 
chapter  2,  for  the  perturbed  Nyquist  plot  not  to  change 
encirclements  is  that  the  norm  of  the  perturbation  £  £ 
remain  less  than  the  norm  of  the  return  difference  matrix  as 
expressed  in  equation  3.2. 

I  I  &  2  CM  ||<  1/1  |  (3+fi}-»)  I!  (3.2) 

cJ  >0  Ibis  condition  will  guarantee  that  the  Locus  of  the 
det  (I+g)  does  not  pass  through  the  -1  point.  If  the  1L  or 
Euclidean  norm  is  assumed  for  this  condition  the  equation 
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3.2  may  te  expressed  in  terns  of  singular  values  as  equation 
3.3. 


cr  <4  G)<  <T  (I  ♦  G)  (3.3) 

This  result  states  that  as  long  as  the  maximum  singular 
value  of  the  perturbation  matrix  A  §  is  below  the  minimum 
norm  value  cf  the  return  difference  matrix  the  system  will 
remain  stable.  The  problem  of  guaranteeing  robustness 
becomes  that  of  finding  the  largest  norm  of  the  perturbation 
quantity,  the  largest  singular  value,  for  which  the  smallest 
norm  or  singular  value  of  the  return  difference  matrix  will 
remain  non  singular. 

The  multiplicative  form  for  a  system  such  as  figure  3.4 
gives  the  similar  norm  equation  in  3.4 

i  I  *<i(jw)  IJ<1/I  I  (I  ♦  (G)-‘)“IJ  (3.4) 

OJ  >0.  which  may  te  expressed  as; 

^  (AG)  <  g  (  !  ♦  G-1  )  (3.5) 

Singular  value  decomposition  software  is  readily  available 
to  determine  how  near  the  matrix  I+G  or  1+  (G)“*  is  to 

singularity. 


IV.  OPTIMIZATION 


The  purpose  of  this  chapter  will  be  to  briefly  describe 
several  of  the  currently  employed  optimization  techniques 
and  the  Automated  Design  Synthesis  (ADS)  program,  [Ref.  10], 
which  employs  these  techniques.  In  general,  optimization 
implies  finding  the  "best"  possible  solution  to  a  problem. 
In  actuality  the  best  solution  found  by  an  optimization 
technique  migh+  really  only  be  a  "better"  solution  to  the 
problem.  The  purpose  of  ADS  and  other  optimization  routines 
is  to  allow  a  rational  search  to  be  conducted  to  find  the 
test  possible  design.  The  techniques  of  numerical  optimiza¬ 
tion  are  used  to  logically  vary  the  various  parameters  that 
affect  the  design  until  a  good  solution  is  found. 

As  an  example  of  an  unconstrained  optimization  problem 
consider  the  following  problem  developed  in  [Ref.  11].  The 
problem  is  to  minimize  the  function 

F  (x)  =10x/  4-20x,  aXj  +  lOx**  +x ,  2-2x,  +5  (4 . 1) 

I  (x)  is  often  called  the  objective  function,  the  cost  func¬ 
tion  or  the  penalty  function.  Since  there  are  no  conditions 
imposed  on  the  design  variables,  x#  and  x  t,  and  no  addi¬ 
tional  limits  imposed  on  the  overall  design,  the  problem  is 
considered  to  be  one  of  unconstrained  minimization.  Figure 
4.1  represents  this  problem  in  the  design  space.  From  the 
figure  it  appears  the  optimum  is  near  the  point  1,1. 
Calculus  may  be  applied  to  determine  the  optimum  exactly. 
Taking  the  derivatives 
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Figure  4. 1  Design  Space  for  Example  problem. 


=40x,  3-40x(  *z*2x,  -2=0 


(4.2) 


9?  (x)/3xl=-20x,*  *20x^0  (4.3) 

and  then  solving  the  set  of  equations  it  is  found  that  x, 
=1.0  and  x  ^  =1 . 0 . 

If  design  conditions  are  imposed  on  the  problem  then  the 
optimization  becomes  one  of  constrained  function  minimiza¬ 
tion.  In  other  words,  while  the  minimum  of  a  function  is 
still  sought,  this  minimum  must  exist  within  the  limits 
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imposed  by  the  design  conditions  or  constraints.  The  minimum 
cf  an  unconstrained  function  will  not  necessarily  he  the 
same  minimum  for  a  constrained  function. 

As  another  example  ,  also  from  [fief.  11]#  a  design  is 
sought  which  gives  the  minimum  weight  of  a  particular 
column.  This  weight  is  expressed  as: 

W=  </«Ah=  prz  Dh  (4.4) 

where  is  the  unit  weight  of  the  material  and  A  the  cross 
sectional  area.  The  stress  in  the  column  is  given  by: 

CT=P/A  =  P//T  Dt  (4.5) 

without  going  into  detail  the  design  will  be  constrained  by 
the  allowable  stress  on  the  structure.  Other  constraints  , 
Euler  buckling  and  shell  buckling,  are  also  of  interest  to 
the  designer  of  this  column.  The  design  problem  is  then 
stated  as: 

minimize  W  =  Eth  (4.6) 

for  the  constraints 

g  (1)  =  °  /&  -  1<  0  (4.7) 


g<2)  =  &AT0  -IS  o  (4.8) 


9  O)  *  “ISO  (4.9) 

with  C>  10e-06  and  t  >10e-0  6. 
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g  (4)  =  t-D<  0 


(4.  10) 


5s  »C Tj  *  refer  to  tiie  allowable  design  stress,  Euler 

buckling  stress  and  the  shell  buckling  stress  respectively. 


Wall  UiictncM  ratio,  t/t0 


Figure  4.2  Design  Space  for  Column  Problem. 


Figure  4.2  is  a  design  space  diagram  for  this  problem.  It 
is  noted  from  the  figure  4.2  that  the  optimum  is  not  unique 


and  can  be  any  value  along  the  arc  noted  as  optimum  in  the 


The  general  optimization  problem  is  then  written  as: 
minimize  F(x)  objective  (4.  11) 


subject  to  g  •  (x)  S  0  j=1,  m  inequality  constraint  (4.12) 

J 

hK(x)  =0  k=1  ,1  equality  constraint  (4.13) 

x^<  xt-<  x?  i=1fn  side  constraint  (4.  1 4 ) 

where  X  =  col  (x1,x2,.. .  ,xn)  is  the  design  variable. 

The  methods  used  to  solve  this  problem  are  usually 
iterative.  After  the  establishment  of  an  initial  set  of 
variables  the  optimizer  will  update. this  initial  value  until 
the  optimum  values  are  found.  Again  , borrowing  from 
[Ref.  11],  the  iterative  technique  may  be  demonstrated  by  a 
simple  example.  Figure  4.3  is  used  to  illustrate  this 
problem.  Given  the  initial  data  set  X°  the  formula 

xt  a  JO  +  **  s  i  (4.  15) 

can  be  used  to  upgrade  this  estimate  of  X.  The  vector  S  is 
the  search  direction  for  the  iteration  and  the  scalar  quan¬ 
tity  is  the  distance  of  the  move  in  the  S  direction. 
Eeginning  at  X°  it  is  desired  to  reduce  the  objective  func¬ 
tion.  The  search  for  values  of  X  that  reduce  the  objective 
function  is  made  in  the  S  direction  which  in  this  example  is 
the  opposite  of  the  gradient  of  the  function  at  point  X°. 
The  choice  cf  S  could  be  arbitrary  as  long  as  it  reduces  the 


r. 


Figure  4.3  Illustrative  Example  for  Iteration. 

value  of  the  objective.  If  S  is  chosen  to  be  opposite  the 

gradient  of  the  objective  then  the  search  would  be  a 

steepest  descent  search.  Once  the  direction  is  chosen  then 
the  scalar  must  be  found  to  minimize  the  objective  along 
this  direction  vector  without  crossing  a  constraint 

boundary.  The  technigue  for  finding  p<  is  to  evaluate  I,  the 
objective  and  the  constraint  functions  for  values  of  using 
computer  code  and  numerical  interpolation  to  estimate  <x  . 
This  value  then  gives  the  minimum  value  for  F  (X)  along 

this  search  vector.  New  that  s<  for  this  direction  is  known 
the  new  equation  becomes: 

x*  =  x*  ♦  **S2  (4.  16) 

and  now  X1  is  used  to  start  a  new  search  in  direction  22 
and  compute  a  new  ^*that  continues  to  reduce  the  objective 
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function.  These  one  dimensional  searches  are  continued  until 
no  more  reduction  in  the  objective  can  be  found.  At  this 
point  no  further  design  improvement  is  possible. 

Optimization  techniques  do  not  always  lead  to  the  abso¬ 
lute  optimum  when  applied  to  problems  of  practical  interest. 
The  reasons  for  this  could  be  numerical  ill-conditioning  of 
the  problem  formulation  or  simply  that  there  are  multiple 
solutions  to  the  problem.  Because  of  these  difficulties  it 
may  be  adviseable  to  choose  several  different  starting 
points  for  the  optimization  process  and  use  engineering 
judgement  as  to  the  design  most  applicable  to  the  problem 
under  analysis. 

Considering  the  uncon  strained  case  first  where  the 
desire  is  to  minimize  the  function  F  (x)  ,  it  is  well  known 
that  F  (x)  will  have  a  minimum  where  the  gradient  of  F(x)  is 
zero.  That  is: 

grad  (F  (x)  )  =0  (4.  17) 

with  the  VF(x)  defined  as: 

grad  (F  (x)  )  =  (  (x)/9xt  r  . . .  ,  ^F  (x)  /Ox„)  •  (4.18) 

Figure  4.4  shows  why  this  is  a  necessary  condition  but  dees 
not  guarantee  a  global  minimum.  The  gradient  of  F  (x)  is 
zero  at  all  four  points  A,B,C,and  3.  However,  only  A  and  D 
are  minima.  A  would  be  the  global  minimum  for  the  function 
as  defined  here.  D  would  be  only  a  relative  minimum.  To 
check  that  the  zero  gradient  corresponds  to  a  minimum  the 
Hessian  matrix, i.e  the  matrix  of  second  partial  derivatives, 
can  be  examined  for  positive  definiteness.  A  positive  defi¬ 
nite  Hessian  ensures  a  relative  minimum.  The  only  way  to 
prove  a  global  minimum  for  the  function  is  to  show  that  the 
Hessian  matrix  is  positive  definite  for  all  design  variables 
X.  A  test  that  is  seldom  possible  to  perform. 


Figure  4.4  Plot  of  Various  Points  With  Zero  Gradient. 

If  the  problem  of  minimization  is  a  constrained  problem 
the  situation  is  different  from  that  discussed  above.  Ihe 
objective  function  gradient  does  not  have  to  be  zero  at  the 
optimum.  Figure  4.5  illustrates  this  case.  Using  figure  4.5 
and  assuming  a  start  at  point  A  it  is  necessary  to  choose  a 
search  direction  S  which  will  reduce  the  objective  function 
while  not  violating  the  constraint  functions.  Any  direction 
that  will  reduce  the  objective  is  said  to  be  useable.  This 
is  seen  to  be  the  half  plane  sector  below  the  tangent  to 
7  (x)  at  point  A.  If  a  half  plane  to  the  right  of  the 
tangent  to  the  active  constraint  at  point  A  is  considered 
then  the  feasible  sector  is  formed.  A  combination  of  the  two 


conditions  gives  the  useable,  feasible  sector  where  the 
search  direction  S  must  be  chosen.  In  mathematical  nomencla¬ 
ture  the  above  argument  may  be  stated  as: 

useable  direction  grad  (F  (X) )  *S<  0  (4.19) 


feasible  direction  grad  (  g^(x))»  S<  0  for  all  j  (4.20) 
for  which  g^  (x)  =  0. 

Feint  B  in  figure  4.5  shows  a  point  where  the  gradient 
of  the  objective  and  constraint  point  in  exactly  the  oppo¬ 
site  direction.  At  a  point  such  as  B  the  only  search  vector 
S  that  meets  requirements  for  useability  and  feasibility  is 
tangent  to  the  constraint  boundary  and  to  a  line  of  constant 
objective  function.  This  condition  is  stated  as: 

grad(F(x)  )+Z>igrad(gj  (x)  )  +  2>^|tgrad(hx(x)  )  =  0  (4.21) 

where  Aj  >  0  and  A  -unrestricted. 

Hith  this  brief  background  in  optimization  complete  the 
Automated  Design  Synthesis  (ADS)  program  will  be  briefly 
discussed.  More  detail  on  the  ADS  routine  can  be  found  in 
£Ref.  10].  This  code  was  developed  as  a  follow-cn  to  the 
successful  CONMIN  code  £Ref.  12]  developed  by  7a nderplaats. 
It  is  designed  as  a  black  box  optimizer  which  allows  the 
user  to  choose  combinations  of  one  dimensional  search,  opti¬ 
mization  algorithm  and  optimization  strategy.  These  will  be 
discussed  later.  For  the  user  with  specific  requirements  the 
code  may  be  tailored  by  parameter  modification  to  meet 
specific  requirements.  For  the  work  done  in  this  thesis  it 
was  assumed  the  user  of  the  code  has  no  detailed  knowledge 
of  optimization  and  will  want  to  use  the  code  in  the 


Figure  4.5  Constrained  Optimization  Ezaaple 


simplest  form.  As  such#  the  ADS  calls  made  from  the  main 
program  use  only  default  parameters  and  first  forward  finite 
difference  gradients.  Should  analytical  gradients  be  avail¬ 
able  the  user  could  use  them  within  the  code  if  desired  with 
no  difficulty.  The  calls  made  by  the  user  to  the  ADS 
routines  specify  several  important  aspects  of  the  problem 
solution.  The  user  may  call  for  an  optimization  strategy  to 
be  used  in  the  routine.  This  is  not  reguired  and  its  use 


depends  on  the  problem.  These  strategies  are  discussed  in 
[Bef.  10].  Two  strategies  used  most  often  in  the  analysis 
done  for  this  thesis  have  been  seguential  unconstrained 
minimization  using  guadratic  exterior  penalty  function  and 
the  augmented  Lagrange  multiplier  (AIM)  method.  Others 
available  include  seguential  linear  programming  and  seguen¬ 
tial  guadratic  programming. 

The  basic  optimizer  is  also  chosen  by  the  user  from  two 
unconstrained  and  three  constrained  optimization  algorithms. 
The  unccnstr aine d  algorithms  are  Fletcher-Reeves  conjugate 
directions,  Davidon-Fletcher-Eowell  (DFP)  variable  metric 
method  and  the  Broydcn-Fletcher-Goldf arb-Shanno  (BFGS)  vari¬ 
able  metric  method.  The  method  of  feasible  directions  and 
robust  feasible  directions  are  available  for  constrained 
minimization. 

The  user  has  available  several  types  of  one  dimensional 
searches  using  Golden  Section  or  polynomial  approximation 
technigues.  The  ADS  code  has  tailored  these  one  dimensional 
search  algorithms  for  the  unconstrained  and  constrained 
cases,  allowing  the  user  to  make  appropriate  choices  for  the 
type  of  problem  to  be  solved.  Figure  4.6  shows  the  basic 
organization  of  the  ACS  program. 

Since  it  is  not  the  purpose  of  this  chapter  to  cover  the 
subject  of  optimization  in  detail  all  the  possible  routines 
in  ADS  will  not  be  discussed.  A  few  of  the  routines  found  to 
work  well  for  the  wcrk  completed  in  this  thesis  will  be 
briefly  described,  however.  There  are  several  methods  of 
optimizing  functions  cf  one  variable  or  one  dimensional 
searches.  For  instance  a  large  number  of  points,  n,  could  be 
chosen  and  the  function  F(x)  evaluated  at  each  point.  The 
point  corresponding  to  the  smallest  value  of  F  (&)  could  then 
be  considered  the  optimum  value  of  the  one  dimensional 
search.  This  method  is  hit  and  miss  and  better  methods  of 
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USER-SUPPLIED  CALLING  PROGRAM 


Figure  4.6  Organization  of  &DS  Program. 

locating  optimum  points  can  be  used.  The  methods  used  in  ADS 
are  Golden  Section  ard  polynomial  interpolation. 

Golden  Section  search  methods  are  easy  to  program  on  the 
digital  computer  and  do  not  require  continuous  derivatives. 
They  have  a  known  convergence  rate  and  are  reliable  for 
poorly  conditioned  problems.  The  major  drawback  of  the 
Golden  Section  routines  is  the  large  number  of  function 
evaluations  require.  The  Golden  Section  algorithm  is  simply 
illustrated  through  the  use  of  figure  4.7  Assume  that  X°  and 
X1  are  known  to  be  bounds  on  the  curve's  minimum  value. 
Also,  the  function  values  of  F  (x)  ,  F°  and  F1  ,are  evaluated 
and  known  at  these  points.  By  picking  two  intermediate 
points  X2  and  X3  where  X2  <  X3  and  evaluating  the  function 


Figure  4.7  Golden  Section  Diagram. 

at  these  points  a  hound  on  the  minimum  may  be  modified. 
Since  F2  for  this  figure  is  larger  than  F3  the  point  X2 
forms  the  new  lower  bound.  The  minimum  is  now  between  X2  and 
X1  .  If  the  function  F*  at  X4  is  determined  and  shown  to  be 
larger  than  F3  then  X4  becomes  the  new  upper  bound  on  the 
minimum.  Ey  repeating  this  procedure  the  bounds  may  be 
narrowed  to  any  desired  tolerance.  The  Golden  Section  rule 
is  applied  to  this  problem  to  reduce  the  bounds  in  the 
quickest  possible  time.  By  appropriately  picking  the  values 
of  X's  at  which  each  function  evaluation  is  made  an  effi¬ 
cient  algorithm  that  uses  the  ratio  proportion  of  the  Golden 
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Section, i.e.  X2/Xl  =  1.62803  is  developed.  By  choosing  a 

value  'T  based  on  this  Golden  Section  rule  where  “2"  =  0.38197 
estimates  for  interior  points  X2  and  X3  can  be  made  as: 


X2  =  (1  -  t  )  x°  *  £-X* 


(4.  22) 


X3  =  ?  x°  ♦  (1  -  't)  XI 


(4.23) 


As  each  new  bound  is  found  the  process  repeats  until  the 
accepted  level  of  convergence  is  reached. 

The  polynomial  interpolation  method  is  accomplished  by 
first  fitting  a  polynomial  curve  to  the  points  about  where 
the  micimum  is  desired  and  then  finding  the  minimum  of  the 
polynomial  function.  For  example  if  the  function  F  (X )  is 
approximated  by  a  guadratic  as: 


F  3  a  a  +a  (  X  +  a  t  X2 


(4.  24) 


Then  the  value  of  X,  X  ,  where  F'  is  zero  can  be  shewn 
to  be  : 


X*  =  -a ,  /2a  z 


(4.25) 


If  a  2  is  positive  then  F  will  be  minimum.  Other  degrees  of 
polynomials  may  be  used  in  similar  fashion. 

Now  that  the  basic  one  dimensional  search  methods  have 
been  reviewed  the  next  step  is  to  examine  the  basic  opti¬ 
mization  routines.  First,  the  unconstrained  case  will  be 
reviewed.  The  optimum  X*is  at  the  point  where 
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grad  (  F  (/) )  =  0 


(4.26) 


Several  zero  order  methods  exist  for  the  purpose  of  finding 
the  minimum  value.  Ihese  include  random  search,  Powell’s 
method  and  Box’s  method.  Since  ADS  does  not  use  these 
methods  they  will  not  be  discussed. 

Automated  Design  Synthesis  makes  use  of  first  order 
methods  which  will  new  be  discussed.  The  steepest  descent 
method  is  best  known,  but  poor  in  performance.  Just  as  the 
name  implies  a  search  direction  is  chosen  opposite  the 
gradient  of  the  objective  function.  At  iteration  0 

S o  =  -  grad  (  F  (X°) )  (4.  27) 

Figure  4.8  shows  this  algorithm  geometrically.  Note  that 
the  method  simply  stair  steps  its  way  down  the  "hill”  to  the 
valley  or  minimum.  ACS  uses  the  Fletcher-Reeves  modification 
to  steepest  descent.  In  this  routine  a  conjugate  direction 
is  chosen  to  improve  the  speed  of  the  search.  Figure  4.9 
shows  how  this  method  tracks  to  the  solution. 

The  variable  metric  methods  listed  earlier  are  usually 
more  powerful  than  Fletcher-Reeves  because  they  store  infor¬ 
mation  that  allows  the  algorithm  to  approximate  the  inverse 
of  the  Hessian  matrix  or  second  derivative.  For  further 
discussion  cf  these  methods  see  [Ref.  11]  or  other  similar 
optimization  texts. 

AES  employs  two  direct  methods  for  constrained  minimiza¬ 
tion.  One  method  is  that  of  feasible  direction  and  the  ether 
is  the  method  of  robest  feasible  direction.  Since  these  two 
methods  were  seldom  employed  within  the  work  presented  in 
this  paper  they  will  not  be  discussed. 

The  methods  chosen  to  handle  the  constrained  minimiza¬ 
tion  problems  formulated  in  this  thesis  are  referred  to  as 
Sequential  Unconstrained  Minimization  Techniques  (S DMT) . 
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Figure  4.8  Steepest  Descent  Algorithm. 

SUMT  methods  are  methods  which  formulate  the  objective 
function  and  the  constraint  functions  into  an  augmented 
objective  function  and  then  solving  the  problem  as  if  it 
were  an  unconstrained  optimization  task.  &  DS  employs 
interior  and  exterior  penalty  function  techniques  as  well  as 
an  Augmented  Lagrange  Multiplier  (ALM)  technique. 

The  exterior  penalty  function  method  is  incorporated  by 
forming  a  penalty  from  the  constraint  equations.  This 
penalty  is  of  the  form  : 


P  (*)  =  21  max  (0,g  .  (x))**^  (h*(x))  2 


(4.  28) 


*2 


Figure  4.9  Fletcher-Beeves  Conjugate  Direction  Algorithm. 

From  this  eguation  it  is  obvious  that  P  (x)  is  zero  if  all 
the  constraints  are  satisfied.  That  is,  if  g  .  (x|  <  0  and  k 
(x)  =  0,  then  all  conditions  of  the  penalty  function  are 
satisfied.  If  an  element  of  the  penalty  function  is  violated 
then  the  penalty  increases  as  the  square  of  the  violated 
constraint.  A  pseudo  or  augmented  objective  function  is 
formulated  where  : 

9  FU)  ♦  r^  P  (x)  (4.  29) 

The  constant  r^  is  a  weighting  factor  for  the  penalty.  It  is 
adjusted  with  ADS  as  the  optimization  proceeds  to  allow  the 
program  to  systematically  converge  to  an  optimum  solution. 
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One  disadvantage  of  this  method  is  that  if  the  optimizer  is 
stopped  short  of  the  optimum  the  design  will  he  in  the 
infeasible  region  and  probably  not  useful. 

Originally  developed  as  a  method  to  solve  equality 
constrained  problems  the  version  of  ALM  in  ADS  has  been 
modified  to  work  with  both  equality  and  inequality 
constraints.  The  statement  is  as  follows: 

minimize  F(x)  (4-30) 


subject 

ho  hK  (x)  =  0  k=1  ,1 

(4.31) 

Next  a  Lagrangian  is  created  such 

t'hat 

L  (x,» 

1  . 

=  F(x)  ♦  ZAk  hA(x) 

(4.32) 

Then  the 

eguation  is  augmented 

with  an  exterior 

penalty 

function  such  that 

A  (x,  h  . 

rp)=  F(X)  ♦51UllhK(x)  + 

rp(h^(x))z  ) 

(4.  33) 

The  power 

of  this  method  is 

that  in  theory 

precise 

constraint 

matching  is  possible 

whereas  in  the 

exterior 

penalty  function  method  it  is  not.  The  full  details  of  the 
method  will  not  be  covered  here,  however,  the  final  form  of 
the  objective  function  will  be  included: 

+  +rp(h  (x))2  )  (4.  34 ) 

j  »  i  Hzi  *►«»  r 

where  =maz  (g  j  (x)  ,  ->i/2rP  } 

the  Lagrange  multipliers  are: 


9 


and  the  update  formulas  for 


(4. 35) 


P 

♦  2rp  (max  (g^  <x)  ,->j/2rp ) )  j=1,m 

C  ^L*  i"1-1  t4-36) 

With  this  brief  review  of  optimization  and  the  ADS  program 
completed  the  next  chapter  will  discuss  the  progcam  develop¬ 
ment  for  this  thesis  research. 
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¥.  OPTIMIZATION  DESIGN  PBQCEDOBE 


Tc  accomplish  the  task  of  designing  a  control  system 
with  acceptable  time  domain  performance  and  robustness 
characteristics  in  a  straight  forward  manner  a  numerical 
optimization  technique  was  chosen  as  the  method  of  iiplemen- 
tation  of  the  design  algorithm.  Using  numerical  procedures 
to  adjust  selected  design  variables,  in  this  case  the  feed- 
hack  and/or  filter  gains,  a  desired  level  of  performance  can 
he  achieved.  This  level  of  performance  is  actually  a  combi¬ 
nation  of  time  domain  performance  and  robustness  or 
frequency  domain  performance.  By  establishing  the  criteria 
for  the  system  performance  in  terms  of  an  optimization 
objective  and  constraint  functions  a  versatile  procedure  can 
he  developed  to  set  the  system  feedback  gains  and  affect  an 
acceptable  design  in  terms  of  performance  and  robustness. 
The  pcle  placement  and  robustness  (POPLAB)  design  procedure 
uses  pole  placement  to  establish  a  designer  selected 
performance  level  and  then  a  minimum  singular  value  level  to 
establish  the  robustness. 

The  pole  placement  portion  of  the  procedure  will  be 
discussed  first.  The  pole  placement  technique  was  chosen 
because  it  was  relatively  easy  to  implement  through  a  numer¬ 
ical  optimization  routine.  By  using  this  numerical  procedure 
it  is  also  simple  to  incorporate  robustness  into  the  proce¬ 
dure  along  with  the  performance  requirements.  A  numerical 
technique  similar  to  one  posed  in  [Ref.  13]  was  chosen  for 
the  pole  placement  algorithm.  An  unconstrained  optimization 
routine  frca  the  IBM  IMSL  library  was  used  for  this  program. 
The  routine,  a  Newton  method,  uses  adjustments  to  the  output 
feedback  gains  to  reduce  the  size  of  an  objective  function. 
This  objective  function  was  expressed  as  a  function  of  the 


pole  location  of  the  system  and  as  the  objective  was  reduced 
the  poles  were  moved  toward  the  desired  locations.  To 
provide  more  versatility  in  the  pole  placement  algorithm  a 
method  that  can  use  constraint  functions  as  well  as  uncon¬ 
strained  optimization  was  chosen  for  this  program.  Ihe 
designer  may  use  either  objective#  constraint  or  a  combina¬ 
tion  of  functions  to  secure  the  desired  pole  locations.  As 
currently  implemented  in  the  program  the  cost  or  objective 
portion  of  the  pole  placement  procedure  is  constructed  as 
equation  5.1 


(5.1) 


OBJ  =  Z  (Av-Vfci  )2  +  <Aj "Ali  )2  (5.1) 

<-/  *  * 

where  =  real  eigenvalue 

=  imaginary  eigenvalue 
\#  =  desired  eigenvalue  location 
=  desired  eigenvalue  location 
The  constraint  formulation  is  a  function  ti at  must  be 
kept  negative  or  the  constraint  is  violated.  It  is  written 
a's  eguation  5.2 


g(j)  =  V  )2  ♦ 


2  -  r 


(5.2) 


where  r  is  a  tolerance  circle  established  as  a  function  of 
pole  placement  position.  Since  the  aim  of  the  optimizer  is 
to  keep  g  negative  any  time  the  ^  function  of  the 
constraint  is  greater  than  r  the  constraint  will  become 
active#  i.e.  violated.  The  optimizer  will  then  attempt  to 
move  the  constraint  to  the  inactive  status  by  adjusting  the 
design  parameters  of  the  system. 

Consideration  of  implementation  of  the  frequency  domain 
or  robustness  porticn  of  the  design  procedure  begins  with 
the  concept  of  MIMO  phase  and  gain  margins.  Several  useful 
theorems  on  singular  value  analysis  of  multiloop  systems  are 


presented  in  [Sef.  8].  One  of  these  theorems  relates  the 
matrix  singular  value  of  the  return  difference  function  to  a 
parameter,^,  and  further  shows  that  as  long  as  the  maximum 
singular  value  of  the  perturbation  function  (L  -  I) 
remains  less  than  this  the  system  remains  stable.  The 
value  of  o(  is  then  related  to  gain  and  phase  margins  of  the 
MIMO  system.  The  relationship  developed  is  given  in  equa¬ 
tions  5.3  and  5.4: 

gain  margin  =  GM  =1/(1  +  *o)  (5.3) 


phase  margin  =  PM  =  icos  ‘1  ( 1-*02/2)  (5.4) 

provided  that  equation  5.5  holds. 

G“(I  +  G)>  Mo 
for  seme  £<„<  1 

These  phase  and  gain  margins  are  guaranteed  in 
simultaneously . 

Universal  gain  and  phase  margin  curves,  [  Bef 
on  the  minimum  singular  values  of  the  return 
matrix  are  developed  from  eguation  5.6. 

(I-i-I)  =max</  (1-  1/}s^)2  +  2/krt(1-cos4(,rt)  (5.6) 

for  all  n  with  k„  >  0.  These  curves  shown  in  figure  5.1 
allow  the  designer  to  pick  a  singular  value  that  corresponds 
to  a  specific  gain  and  phase  margin  for  a  given  system.  In 
addition  to  the  universal  gain  and  phase  plot  [Bef.  15] 
developes  an  optimizer  solution  for  formulating  a  robust 
controller  using  the  CCNMIN  algorithm  [Ref.  12]. 


(5.5) 

every  loop 

.  14],  based 
dif f ere nee 


Figure  5.1  Universal  Gain  and  Phase  Singular  Value  Plot. 

Since  the  universal  curve  in  figure  5.1  provides  a 
convenient  method  of  specifying  gain  and  phase  margins  in 
terms  of  singular  values  the  robustness  portion  of  the  pole 
placement  and  robustness  design  procedure  uses  the  minimum 
singular  value  level  of  the  return  difference  matrix  to 
determine  the  robustness.  The  minimum  singular  value  level 
is  formulated  as  an  objective  or  constraint  function  in 
eguation  5.7 

J=£.<aax  (0,  -  <T  (>,p)})  2  (5.7) 

J 

The  optimization  procedure  may  be  used  to  change  feedback 
gains  until  the  minimum  singular  value  is  raised  above  this 

desired  design  le-'el.  Although  the  same  formu _ tion  can  be 

used  as  a  negative  constraint  function  it  has  not  been 
implemented  as  such  within  this  program.  There  are  numerous 
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ways  the  singular  value  formulation  could  be  i mplemented 
within  the  program  by  changes  of  the  code  if  design  require¬ 
ments  forced  such  changes. 

The  pcle  placement  and  robustness  design  program  is 
based  on  the  ADS  code  discussed  in  Chapter  4  to  implement 
the  design  variable  selection  procedures.  The  pole  placement 
and  robustness  program  consists  o£  two  separate  programs. 
The  first  program  is  used  to  provide  designs  for  state  or 
output  feedback  problems  while  the  second  program  is  used 
for  observer  or  filter  designs. 

For  the  state  or  output  feedback  design  program  the  user 
must  input  the  plant  matrices  A,B,C  and  initial  starting 
values  for  the  feedback  matrix  F.  The  matrices  correspond  to 
the  following  linear  differential  system: 

x  =  Ax  +  Bu  (5.8) 


y 


=  Cx 
— 


(5-9) 


u  =  -Fx  (5.  10) 

A  feed-forward  matrix  has  not  been  considered  in  the  current 
program.  A  feed  forward  matrix  could  be  added  to  the  proce¬ 
dure  if  required  for  specific  design  cases. 

As  the  design  program  is  currently  coded  the  user  may 
run  output  feedback  or  state  feedback  design  by  specifying 
the  C  matrix  as  the  diagonal(I)  matrix  for  state  feedback. 
The  program  relies  on  initial  starting  values  of  the  feed¬ 
back  gains,  F.  As  discussed  in  Chapter  4  there  is  no 
guarantee  that  the  optimum  found  by  the  procedure  each  time 
is  the  global  optimum  or  that  the  procedure  will  always 


converge  to  an  acceptable  solution.  The  ability  tc  select 
acceptable  starting  values  for  the  feedback  gains  will  make 
the  procedure  more  efficient  in  operation.  As  currently 
employed,  the  program  is  used  to  obtain  pole  placement  and 
robustness  for  a  given  set  of  starting  gains  and  a  selected 
optimization  routine  from  the  ADS  program.  If  the  optimizer 
is  net  able  to  meet  the  desired  design  goals  on  this  program 
run  two  options  are  available.  First,  change  to  a  different 
optimization  routine  from  the  list  of  available  ADS  routines 
and  rerun  the  problem.  This  was  usually  successful  in 
improving  the  design.  Second,  the  designer  uses  a  new  set  of 
starting  values  for  the  feedback  gains  and  repeats  the 
design  procedure.  Both  options  might  be  used  on  particularly 
difficult  cases. 

The  pole  placement  and  robustness  design  procedure  has 
consistently  been  able  to  find  improved  designs;  however  the 
program  does  not  always  yield  acceptable  designs .  Certain 
problems  reguire  changes  in  the  optimizer  routine  and  modi¬ 
fication  in  the  initial  feedback  gain  starting  values  in 
order  to  obtain  acceptable  designs.  Using  the  IBM  3033  time 
share  system  the  pole  placement  and  robustness  routine 
requires  about  ten  CPU  seconds  to  work  a  second  order 
problem  and  on  the  order  of  15  to  60  CPU  seconds  to  run  a 
forth  order  problem.  The  actual  amount  of  time  varies  with 
optimization  requirements  and  time  share  utilization. 

The  observer  robustness  design  program  requires  two 
passes  of  the  ADS  program.  In  the  first  pass  the  feedback 
gains,  F,  of  the  controller  are  computed  to  obtain  the 
desired  pole  locations.  The  second  pass  of  the  ADS  routine 
is  used  to  adjust  the  observer  gains  to  recover  the  system 
robustness.  The  two  pass  procedure  was  chosen  because  it 
allows  a  smaller  number  of  design  variables  at  each  stage  of 
the  optimization  and  much  more  efficient  computer  usage. 
Figure  5.2  shows  how  the  observer  is  implemented.  This 
diagram  is  algebraicly  stated  as; 


X  =  |a  -BK*j  fxl  ♦  [b]  r  (5.  11) 

$  FC  AC-BK  J  [x(  (b[ 

Ae  =  A  -  FC 

where  x  is  the  state,  x  is  the  estimator  variable,  F  is  the 
feedback  gain  and  K  is  the  observer  gain.  The  design  cf  the 
feedback  gains  and  the  observer  gains  are  accomplished  as 
separate  guantities  in  keeping  with  the  separation  prin- 


Figure  5.2  Observer  Implementation. 

ciple.  In  using  the  pole  placement  and  robustness  design 
procedure  for  the  observer  system,  initial  values  of  the  F 
and  K  matrices  must  be  input.  The  same  or  different  opti¬ 
mization  techniques  from  ADS  may  be  employed. 

The  routines  contained  in  this  thesis  have  been  based  on 
the  input  additive  singular  value  level.  The  pole  placement 
and  robustness  design  algorithm  computes  input  additive. 
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output  additive,  input  multiplicative,  and  output  multipli¬ 
cative  singular  values.  Any  of  these  singular  values  can  be 
incorporated  into  tie  objective  or  constraint  formulations 
tut  have  net  been  for  this  version  of  tne  program. 

In  summary,  the  pole  placement  and  robustness  design 
procedure  is  a  straight  forward  numerical  optimization 
procedure  fer  the  practical  application  modern  MIMO  system 
analysis.  The  new  aspects  of  the  procedure  are  the  implemen¬ 
tation  of  both  pole  placement  and  robustness  criteria  within 
the  same  design  program.  The  versatility  of  the  pole  place¬ 
ment  and  robustness  design  is  obtained  by  incorporating  a 
state  of  the  art  optimizer  routine  ADS,  with  currently 
available  singular  value  computation  routines.  Using  the 
optimizer  format  for  the  pole  placement  and  robustness 
design  gives  the  designer  the  ability  to  modify  variables 
directly  that  affect  both  the  time  domain  or  performance  of 
the  system  and  the  frequency  domain  or  robustness  of  the 
system.  The  numerical  optimization  incorporated  into  the 
pole  placement  and  robustness  design  is  flexible  enough  to 
incorporate  the  double  pass  design  technique  using  the  sepa¬ 
ration  of  the  feedback  and  filter  gains  to  obtain  robustness 
recovery  for  observer  based  controllers. 


VI.  IHTRO DOCTOR!  PROBLEM 


The  purpose  of  this  introductory  problem  is  to  review 
some  results  of  modern  multivariable  robustness  theory  using 
a  simple  problem.  This  same  simple  system  is  then  used  as  a 
test  problem  for  the  optimization  technique  developed  for 
this  thesis.  The  problem  provides  excellent  insight  into  the 
cross-coupling  problem  and  demonstrates  how  effectively  the 
pole  placement  and  icbustness  design  procedure  can  be  in 
dealing  with  the  cross  perturbation  terms. 

The  problem  chosen  for  this  introductory  analysis  comes 
from  £Ref.  9].  Figure  6.1  is  a  diagram  of  this  basic 
system.  In  this  problem  a  simple  plant  is  specified  by  the 
following  linear  system: 


x  1 
x2 
where 
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yl  =  xl  (6.2) 

y  2  =  x2 

A  feedback  compensation  of  the  form  of  equation  6.3  was 
assumed. 
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Which  gives  a  closed-loop  system,  equation  6.4 
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stem  has  eigenvalues  at  -2 


Using  equation  6.5 


,-2 


and  is 


(6.4) 

therefore 
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G=  C(sl-l)-*fi 

the  transfer  matrix  may  be  written  as  equation  6.6. 


(6.5) 


G  =  fl/s+1  b  ,c/s*l] 

[  0  1/s+1  J 

which  gives  the  return  difference  matrix# 


(6.6) 


equation  6.7. 


I  +  G(s)  =  s+2/s+l 
0 


bJt/s+  1 
S+2/S+  1 


(6.7) 


♦  y, 


*y2 


Figure  6.1  Basic  Hulti-input  Multi-output  System. 

The  problem  shows  the  inadequacies  of  classical  methods 
in  establishing  the  robustness  of  the  system.  A  brief  review 
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cf  the  results  will  te  presented.  Using  the  return  differ¬ 
ence  matrix  (1+2) ,  the  determinant  may  be  written  as 
equation  6.8. 


det  (!•*£)  -1=2s*3/  (s*1)  2 


(6.8) 


Therefore,  the  multivariable  Nyguist  diagram  will  be  as 


-1  0 


Figure  6.2  Hultivaria ble  Nyguist  for  2s+3/(s+1)2. 


shown  in  figure  6.2.  The  diagram  does  not  encircle  the 
(-1,0)  point  and  is  indicative  of  a  closed-loop  stable 
system.  Considered  as  a  SISG  system  one  has  a  gain  margin  of 
-1/3  to  co  and  a  phase  margin  of  ±106  degrees.  Under  this 
criteria  one  can  conclude  that  this  is  a  good  design.  This 
will  be  shewn  later  net  to  be  the  case. 

flulti-input  systems  are  often  designed  a  loop  at  a  time. 
Applying  loop  design  to  the  system  the  transfer  function  for 
either  lcop  taken  separately  becomes  equation  6.9 


G=  1/S  +  1 


(6.9) 


Figure  6.4  shows  the  Nyguist  diagram  of  this  problem.  Using 
this  Nyguist  diagram  the  system  is  indicated  to  he  stable 
and  have  phase  and  gain  margins  of  and 

EM=(±180).  This  analysis  does  not  show  the  true  nature  of 
the  robustness  of  the  system.  Since  the  factor  b  is  not  a 
parameter  in  either  of  the  Nyguist  curves  it  plays  no  part 
in  the  stability  determination  using  these  diagrams. 

Using  the  criteria  of  singular  values  discussed  in 
Chapter  3  a  measure  of  the  nearness  to  instability  for  this 
problem  may  be  obtained  by  plotting  the  minimum  singular 
value  of  (I+G)  ,  5* .  (For  numerical  calculations  a  value  of 
=  50  is  assumed).  Figure  6.3  shows  the  plot  of  this 
value  vs.  freguency.  This  gives  a  minimum  singular  value  of 
about  -23  db  or  0.07  1  which  corresponds  to  a  gain  margin  of 
about  0.S3  to  1.08  and  a  phase  margin  of  ±4.1  degrees.  These 
phase  and  gain  margins  are  quite  small  and  are  evident  in 
the  cross-feed  perturbation  problem. 

A  perturbed  system  as  shown  in  figure  5.5  can  be 
produced  which  will  lead  to  stability  problems  with  very 
small  values  of  perturbation.  The  closed-loop  system  if 
perturbed  by  a  small  perturbation,  5/bir  ,  where  b  lt  is  a 
large  number,  will  have  as  a  characteristic  eguation,  equa¬ 
tion  €.10 

(si- A)  =  s2  ♦  4s  +5  (6.10) 

with  the  eigenvalues  of  s=-2±  There  is  one  positive 

root  in  this  solution  and  the  system  is  unstable. 

To  determine  the  nature  of  the  robustness  of  this  system 
the  return  difference  matrix  must  be  considered.  If  the 
return  difference  matrix  of  the  transfer  function  (I+G(j«)) 
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Frequency  (rad/sec) 

Figure  6.3  flinimum  Singular  Value  Plot,  Example  Problem. 

is  nearly  singular  at  some  f  reguency  CO, ,  then  the  multivari¬ 
able  system  will  not  be  robust  with  regard  to  any  modelling 
errors  within  the  system.  This  is  because  any  small  change 
in  G  (j«)  can  then  take  I+jG(jw)  singular  and  the  det  (l*£) 
becomes  zero,  thus  changing  the  encirclements  of  the  Nyguist 
stability  point  and  indicating  a  system  instability.  Using 
equation  6.1  the  pole  placement  and  robustness  design  method 
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Figure  6.4  Nyguist  Diagram  of  1/s*1. 


can  be  demonstrated  .  The  first  step  is  to  establish  a  base¬ 
line  for  the  design.  The  state  unity  feedback  model  in  this 
simple  problem  gives  pole  locations  of  -2  and  -2.  This  set 
of  eigenvalues  -2,  -2  were  chosen  as  the  baseline  for  the 
system.  Since  a  second  order  system  requires  only  two  feed¬ 
back  gains  to  place  the  poles,  the  diagonal  feedback  gains 
were  chosen  for  pole  placement  purposes.  The  pole  placement 
and  robustness  design  program  was  then  used  to  place  the 
poles  and  in  the  process  obtained  the  gains  required  to  do 
this  (1,1).  A  plot  cf  some  of  the  singular  value  criteria 
obtained  is  shown  in  figure  6.6.  Pole-zero  plots  of  the 


Figure  6.5  Perturbed  Systea. 

closed-  loop  transfer  functions  of  the  closed-loop  transfer 
matrix  are  shown  in  figure  6.7.  In  this  case  the  pcles  and 
a  zero  are  clustered  about  the  -2  point  and  in  the  input  two 
to  output  one  channel  a  zero  is  located  at  the  -1  point  on 
the  pole-zero  diagraa.  This  point  closely  corresponds  to  the 
minimum  singular  value  frequency.  The  only  significant 
observations  are  the  ability  of  the  pole  placement  and 
robustness  routine  to  place  the  poles  and  the  relatively 
poor  singular  values  indicative  of  low  robustness. 

Since  there  is  a  requirement  for  two  feedback  parameters 
to  set  the  poles  of  a  second  order  systea  there  should  be  no 
additional  freedom  in  design  to  account  for  robustness. 
Case  two  was  a  run  to  demonstrate  this  fact.  Again  allowing 
only  two  design  variables  for  the  problem  the  pole  placement 
and  robustness  design  program  was  run  but  with  an  objective 
function  formulated  to  adjust  the  singular  value  above  a 
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SINGULAR  VALUE 


Figure  6.6  Singular  Value  Plot  for  Simple  Problem. 
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base  level  of  0.6.  This  level  would  correspond,  using  the 
universal  gain  margin  chart,  to  a  gain  margin  of  -4  db  to  3 
db  and  a  phase  margin  of  ±  35  degrees.  With  only  the  two 
design  variables  to  work  with  the  pole  placement  and  robust¬ 
ness  design  program  was  unable  to  place  the  poles  and  adjust 
the  singular  value  level  to  the  required  value.  The  poles 
were  placed  at  -2,  -2  but  the  singular  value  minimum  was 
still  cn  the  order  of  -23  db  or  0.07.  The  pole-zero  plot 
remained  almost  unchanged.  It  is  clear  that  additional 
degrees  of  freedom  for  the  pole  placement  and  robustness 
design  program  must  be  opened  if  robustness  is  to  be 
accounted  for. 

Case  three  was  then  run  on  the  pole  placement  and 
robustness  program  by  adding  on  additional  degree  of 
freedom.  This  case  used  the  f  gain  as  the  additional 
design  variable.  A  good  choice  as  will  be  seen.  Allowing 
the  optimizer  routine  the  extra  freedom  to  adjust  the  addi¬ 
tional  feedback  gain  term  an  excellent  design  was  fcund. 
The  singular  value  minimum  became  0.  385  as  shown  in  figure 
6.3.  Using  the  universal  gain  margin  chart  this  corresponds 
to  gain  and  phase  margins  cf  -6  db  -  18  db  and  ±52  degrees 
respectively.  This  is  a  considerable  improvement  ever  the 
original  design.  The  factor  that  changed  the  design  was  the 
upper  diagonal  feedback  term  which  provides  a  cancelling 
factor  for  the  cross-coupling  term  b =  50.  This  can  be 
seen  by  looking  at  the  system  matrix  equation  6.1  1. 

=  A  -  BF  (6.11) 
where  equation  6.12  gives  A 

~  ^  n  “b  it  (6.12) 


» 
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Figure  6.8  Singular  Value  Plot  Case  Three. 
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If  the  upper  diagonal  term  of  the  matrix  can  be  driven  to 
zero  the  system  will  be  decoupled  into  a  diagonal  system. 
To  place  the  poles  cf  this  system  at  -2,-2  the  diagonal 
feedback  gains  must  be  one.  The  optimizer  feedback  gains 
for  this  case  were  fn  =1.00  006,  ^^=0. 99998  and,  most  impor¬ 
tantly,  f/it  =-5  1. 8340  1.  Using  these  values  in  the  A  matrix 
A  =  -2.00006  1.83501 

0.0  -  1.99998 

The  value  in  the  upper  right  position  in  the  system  matrix 
has  teen  lowered  considerably  from  the  value  of  near  50  that 
appears  in  this  position  in  the  low  singular  value  cases, 
lowering  this  system  gain  value  decreases  the  cross-coupling 
perturbation  effects  cn  the  system.  This  can  be  more  graphi¬ 
cally  demonstrated  by  figures  6.9  and  6.10.  In  figure  6.9 
the  transfer  function  shows  a  high  gain  of  approximately  35 
db  and  a  bandwidth  of  50  rad/sec.  In  figure  6.10  this  gain 
has  been  reduced  to  6.0  db  with  a  bandwidth  of  almost  2 
rad/sec.  Two  things  are  indicated  by  the  figures,  one,  the 
open-loop  Ecde  plot  of  the  transfer  function  of  the  cross- 
coupled  channel  can  be  used  to  indicate  the  robustness 
problem  as  evidenced  by  the  high  gain  and  bandwidth  relative 
to  the  other  feedback  gains  and,  two,  the  mechanism  used  by 
the  pole  placement  and  robustness  design  procedure  to 
recover  robustness  is  to  reduce  the  relative  gain  and  asso¬ 
ciated  bandwidth  within  the  affected  channel.  Figures  6.11 
and  6.12  which  are  for  the  input  one  to  output  one  channel 
are  approximately  tie  same  as  are  the  Bode  plots  for  the 
input  two  to  output  two  channel  which  are  not  shown.  These 
figures  indicate  that  no  problem  exists  in  the  diagonal  or 
direct  coupling  terms.  The  pole-zero  diagram  for  optimiza¬ 
tion  run  case  three  is  included  in  figure  6.13.  The  only 
significant  change  in  this  plot  as  compared  to  figure  6.7  is 
the  movement  of  the  zero  in  the  off  diagonal  pole -zero  plot. 
The  zero  is  seen  to  shift  to  one  of  the  pole  locations 
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Figure  6.10  Transfer  Function  2-1  Optimized  (  3  Var). 
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during  the  cptim izaticn.  This  zero  shift  has  the  effect  of 
smoothing  the  frequency  response  curve  in  the  vicinity  of 
the  frequency  of  the  minimum  singular  value  providing  a  more 
uniform  gain  distribution  around  this  point.  Table  1  shows 
comparison  results  for  the  feedback  gains  for  saveral  cases 
of  this  basic  problem. 

The  logical  extension  of  the  problem  to  the  general  case 
is  to  allow  all  four  feedback  gains  to  become  design  vari- 

T  ABLE  1 

Comparative  Results  Simple  Problem 


Feedback  Gains 

fll 

f  12 

f  21 

f  22 

1.00006 

-51.83401 

0.0 

0.9998 

0.25254 

0.0 

-0.01071 

1.52332 

.99960 

-50.0 

0.00025 

0.99038 

1.11936 

-60.4718 

0.00280 

0.75478 

ables.  The  pole  placement  and  robustness  routine  can  then 
use  full  freedom  in  choosing  all  four  of  these  feedback 
gains  to  compensate  for  any  cross-coupling  effects  within 
the  system.  Based  on  the  previous  analysis  the  two  diagonal 
gains  would  be  expected  to  approach  1  while  the  upper  off 
diagonal  gain  moves  tc  -50  and  the  lower  one  moves  to  0.  The 
flexibility  of  the  ACS  program  was  required  for  the  four 
design  variable  study.  Several  runs  were  made  with  various 
combinations  of  starting  parameters  for  the  feedback  gains 
and  optimizer  routines  before  a  good  design  for  the  case 


incorporating  all  four  feedback  gains  as  variable  was 
obtained.  While  this  points  out  one  of  the  limitations  of 
optimization  design  routines  the  program  was  able  to  develop 
an  improved  design  ever  the  baseline  while  employing  four 
design  variables.  In  the  design  produced  for  this  case  the 
singular  value  level  was  placed  above  0.6  for  a  gain  margin 
of  -4  db  to  9  db  and  phase  margin  of  ±  35  degrees.  In  this 
formulation  the  optimizer  was  able  to  place  the  poles  and 
meet  the  design  singular  value  level.  The  feedback  gains 
produced  by  the  optimizer  were: 

1.1  1936  -60.4718 

C.0028  0.75478 

which  are  approaching  the  analytic  design  gains.  After 
obtaining  the  feedback  gains  the  OPTSYS  program  was  used  to 
obtain  the  necessary  data  to  do  a  closed-loop  poLe-zerc  map. 
This  plot  is  shown  in  figure  6.14  These  plots  indicate  a 
similar  pole-zero  location  to  that  found  in  the  previous 
three  design  variable  problem.  Again  the  gain  in  the 
affected  channel  has  teen  reduced  to  compensate  for  the  high 
cross-coupling  perturbation  within  the  system. 

The  design  studies  presented  to  this  point  have  been 
based  on  breaking  the  system  loop  at  the  input  as  shown  in 
figure  6.15.  In  multivariable  theory  the  location  of  the 
break  in  the  loop  changes  the  return  difference  for  the 
system  and  the  transfer  function  formulation.  In  the  figure 
number  1  depicts  a  system  with  an  input  loop  break  point 
while  number  2  is  an  output  loop  break  point  for  output 
return  difference  determination.  The  return  difference  func¬ 
tion  for  the  point  1  is  written  as  I  +  Fg  while  the  return 
difference  fer  point  2  is  I+GF.  The  baseline  system  not  only 
has  low  singular  values  of  the  input  return  difference 
matrix,  the  lowest  being  0.0706,  but  also  has  low  singular 
values  for  the  output  return  difference  matrix  of  point  2  of 
about  the  same  order.  To  demonstrate  the  versatility  of  the 
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Figure  6.15  System  Block  Diagram. 

pole  placement  and  robustness  routine  a  design  was  made  in 
which  the  output  singular  values  were  specified  instead  of 
the  input  values.  The  pole  placement  and  robustness  routine 
produced  a  design  with  gains  that  provided  a  highly  robust 
system  on  the  output  side  with  gain  and  phase  margin  cf  -6 
db  to  and  phase  margin  of  ±60  degrees.  The  design  was  not 
robust  on  the  input  side.  Thus,  designing  for  robustness  at 
one  point  in  the  system  does  not  necessarily  give  robustness 
at  all  points  within  the  system. 

Cne  final  case  that  should  be  discussed  is  that  of 
setting  both  input  and  output  robustness  criteria  at  the 
same  time.  Excellent  results  were  obtained  for  this  case. 
The  design  routine  placed  the  poles  at  -1.97±.009j  with 
feedback  gains  of 

0.85793  -45.  91757 

0.00425  0.87137 

The  input  singular  values  were  raised  to  a  level  of  0.74984 
or  -4.5  db  to  12  db  gain  margin  and  ±43  degrees  phase 
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margin.  The  output  values  were  raised  to  above  0.822  which 
corresponds  to  -5  to  15  db  gain  margin  and  ±49  degrees  phase 
margin. 

Tc  summarize,  it  can  be  stated  that  the  robustness 
problem  for  this  system  exists  in  the  upper  cro ss- coupling 
channel  (input  two, output  one).  The  lack  of  robustness  can 
be  discovered  in  two  ways.  The  first  method  is  to  plot  the 
open-loop  Bede  plots  cf  each  element  of  the  transfer  matrix 
and  lock  for  extremely  high  gains  and  bandwidths  relative  to 
the  ether  transfer  functions.  The  second  method  examines 
the  singular  values  of  the  return  difference  matrix  for 
magnitude.  Low  singular  values  correspond  to  low  robust¬ 
ness.  The  pole  placement  and  robustness  design  routine  can 
increase  robustness  by  modifying  feedback  gains  tc  reduce 
the  effect  of  cross-coupling  within  the  system.  Observing 
the  gain  modification  made  by  the  pole  placement  and  robust¬ 
ness  routine  the  critical  channel  within  the  system  that 
affects  the  robustness  may  be  determined  from  the  Bode 
plots.  The  pole  placement  and  robustness  routine  feedback 
gain  changes  also  cause  zero  shifts  during  the  robustness 
recovery.  The  gain  on  the  open  loop  Bode  plot  for  the 
affected  crcss-coupling  channel  is  adjusted  and  the  closed- 
loop  zeros,  as  seen  on  the  pole-zero  diagram,  are  shifted. 
This  zero  shift  is  in  a  direction  which  will  combine  with 
system  poles  to  smooth  the  freguency  response  diagram  in  the 
vicinity  of  the  minimum  singular  value. 
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VII.  A  H EIICOPT EB  STABILITY  PROBLEM 

This  chapter  will  deal  with  a  more  practical  application 
cf  the  numerical  optimization  program.  In  this  problem  the 
combined  pole  placement,  robustness  design  procei ure  will  be 
applied  to  the  linear  lateral  dynamic  channels  of  a  CH-47 
helicopter.  The  model  is  a  highly  coupled  two-input  two- 
cutput  system  that  has  been  studied  for  its  basic  robustness 
characteristics  [Ref.  16].  The  usual  procedure  for  design 
of  highly  coupled  systems  is  to  obtain  a  diagonally  dominant 
closed-loop  system.  This  diagonal  system  will  be  stable  but 
not  rcbust  to  cross-feed  parameters.  Sandell,et  al,  produced 
three  designs.  Two  cf  these  designs,  while  meeting  basic 
performance  criteria,  had  poor  robustness.  The  third  design 
was  a  relatively  gocd  design.  The  numerical  optimization 
technigue  developed  in  the  thesis  was  applied  to  the  two 
poor  designs  and  shown  to  provide  substantial  improvement  in- 
robustness. 

The  systems  were  designed  to  satisfy  specifications  to 
step  input  response  and  stability  margins  as  stated  in  mili¬ 
tary  specifications.  Specific  design  parameters  for  each  of 
the  three  designs  presented  were  not  available.  The  designs 
were  all  considered  to  meet  the  performance  specification 
criteria  and  stability  margin  requirements.  It  was  shown 
that  two  of  these  designs  were  extremely  sensitive  tc  model 
errors.  The  classic  Nyguist  techniques  did  not  predict  this 
sensitivity.  The  singular  value  analysis  did  indicate  sensi¬ 
tivity  problems. 

The  model  is  that  of  a  CH-473  helicopter  lateral  dynamic 
system  in  hover.  The  dynamic  model  of  the  system  is 
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x  =  A  x  +  B  u 
—  /*  “  **  * 


(7.1) 


X  =(v,p,r,y) 


(7.2) 


'4  =  (  &£  t  Sc  ) 


(7.3) 


where 


2.27 

0.01 


-1.42  -0.15 


-0.7 


0.04  -0.05 


and  where 


0.  12 


0.95 


-0.  07 

-0.  5 
0.  11 


31.99 


0.04  -8.37 

.34  0.  0  2 


with  full  state  available  for  feedback.  Table  3  is  a  summary 
of  parameters.  The  system  is  not  open-loop  stable. 

Three  control  laws  are  formulated  to  satisfy  the  desired 
performance  specifications.  Equation  7.4  is  the  basic 
control  law. 


U  =  -F  -X  ♦  h*. 

—  Ml  l  ' 

i~  1,  2  or  3 

where  the  following  values  of  ?  and  h  were  are; 


(7.4) 


f,  - 


£l  = 


-1.72  -23.5  70.6  595. 

0. 024  -2.71  0.368  -7.99 

0.  198  154.0  13.  3  142. 

.-0.01  -1.  592  -0.  139  -1.4 
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TABLE  2 

CH-46  Helicopter  Parameter  Definitions 


Variable 

Units 

1 

Description 

V 

ft/sec 

Vehicle  body  y-axis  earth  relative  velocity  component 

P 

rad/sec 

Roll  rate 

r 

rad/ sec 

Yaw  rate 

4 

rad  | 

Roll  attitude  angle 

63 

in 

Yaw  rate  rotor  deflection  control 

6c 

in 

Roll  rate  rotor  deflection  control 

The  **Pc  is  a  step  input  command  that  must  track.  Figure 
7.1  is  a  diagram  of  the  control  structure  for  the  feedback 
control  laws.  Figure  7.2  is  a  detailed  layout  of  the 
system.  All  three  of  the  designs  have  negative  real  eigen¬ 
values  and  provide  stable  overdamped  responses.  All  three  of 
the  above  designs  meet  the  desired  design  specifications. 

The  loop  Bode  plots,  [fief.  16],  indicate  all  three 
designs  to  have  good  stability  margins  when  considered  a 
loop  at  a  time.  Since  these  designs  were  all  full  state  LQ 
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Feedback 

Controller 

Gain 


designs  with  a  diagonal  control  weighting  matrix  they  should 
posses  at  least  -6  db  to  infinite  gain  margin  ani  60  degrees 
of  phase  margin.  lihen  the  singular  values  of  the  three 
designs  are  computed  a  robustness  problem  is  indicated  by 
low  singular  values  cf  two  of  the  designs.  Figure  7.3  pres¬ 
ents  singular  value  plots  of  all  three  designs.  This  plot 
shows  that  designs  1  and  2  both  have  very  low  minimum 
singular  values  for  the  return  difference  matrix.  Design  1 
goes  as  low  as  -20  db  near  10  rad/sec  in  freguency  while 
design  2  is  down  to  -34  db  at  freguencies  up  to  1000 
rad/sec.  Design  3  is  a  good  design  with  singular  values  that 
remain  above  one  throughout  the  freguency  range  of  interest. 
Using  the  universal  gain  and  phase  diagram  as  discussed  in 
Chapter  6  this  eguates  to  a  gain  margin  of  about  -6  db  to 
infinity  and  a  phase  margin  of  60  degrees  as  expected. 

The  system  stability  in  design  1  may  be  affected  by 
perturbations  occurring  in  the  actuators  as  shown  in  figure 
7.4.  if  the  output  axis  coupling  from  Sc  spills  into  S3  in 
the  freguency  range  from  0.5  to  50.0  rad/sec  with  a  magni¬ 
tude  cf  0.12  and  a  phase  of  60  degrees  then  the  system  can 
become  unstable.  This  could  be  caused  by  nonlinear  terms, 
worn  parts,  or  system  saturation. 

The  second  design  as  shown  in  figure  7.5  may  have  a 
stability  derivative  variation  between  Sa  and  p.  If  this 
derivative  varies  from  about  0.04  to  -0.96  the  system  can 
become  unstable.  The  center  of  gravity  location,  trim  of  the 
aircraft  and  rotor  coupling  can  all  affect  this  derivative. 

The  pole  placement  and  robustness  design  technigue 
developed  for  this  thesis  research  was  applied  to  designs  1 
and  2  to  obtain  an  improved  robustness  for  these  designs. 
Design  1  will  be  considered  first.  In  this  design  a 
cross-feed  perturbation  through  the  actuator  can  produce 
instability  as  indicated  by  the  low  singular  values  at  lower 
freguencies.  To  study  this  problem  with  the  pole  placement 
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and  robustness  design  method  it  is  assumed  that  the  eigenva¬ 
lues  (pole  locations)  of  the  system  as  developed  in 
[Ref.  16].  are  the  required  poles  for  the  performance 
criteria.  Once  the  pole  locations  are  set,  robustness 
criteria  must  be  selected.  From  the  universal  gain  and  phase 
margin  curve  discussed  earlier  two  choices  of  singular  value 
levels  were  made  for  this  problem.  The  first  singular  value 
level  chosen  was  0.6.  This  corresponds  to  a  gain  margin  of 
-4.0  db  to  8  db  and  a  phase  margin  of  about  35  degrees.  The 
second  value  chosen  was  1.0.  With  corresponding  gain  margin 
of  -6  db  to  oo  and  a  phase  margin  of  60  degrees  which  are 
the  characteristics  of  a  LQ  regulator  design  with  diagonal 
weighting  matrices. 

For  singular  value  level  0.6  the  pole  placement  and 
robustness  design  routine  places  the  poles  as  shown  in  table 
2  The  slight  differences  in  these  pole  locations  appear  to 
have  insignificant  effect  on  the  performance  as  shown  in  the 
response  curves  of  the  system.  The  feedback  gain  adjustment 
moves  the  minimum  singular  value  from  about  0.11  with  very 
poor  phase  and  gain  margins  to  a  level  of  0.66  .  This  is 
above  the  desired  values  of  gain  and  pnase.  The  improvement 
in  robustness  came  from  modifying  the  feedback  gains  in 
channel  Ss  .  3y  greatly  reducing  the  gains  in  channel  S&  the 
optimizer  minimizes  the  influence  of  the  cross-coupling  from 
channel  Sc .  In  this  way  a  much  larger  spill  over  of  channel 
Sc  may  be  tolerated  through  the  actuator  before  the  system 
will  become  unstable. 

The  primary  mechanism  of  robustness  improvement  in  this 
problem  was  a  reduction  in  the  gain  levels  of  the  affected 
channel.  The  feedback  gains  presented  in  table  4  show  the 
modification  of  these  gains  from  those  utilized  in 
[Ref.  16].  While  all  the  gains  are  modified  the  gain 
undergoes  a  significantly  larger  change  than  the  ether 
design  1  case  1  feedback  gains.  Loojting  at  tti  e  open-loop 


TABLE  3 

Design  One  Pole  Placement 


Pole 

Desired  Location 

Actual  Location 

1 

-24.7977 

-24.7893 

2 

-11.3635 

-12.0083 

3 

-10.3288 

-10.7752 

4 

-  2.1005 

-  2.1181 

transfer  functions  cf  this  optimized  problem  and  comparing 
them  with  a  non- optimized  set  of  transfer  functions  for  $S 
to  as  shown  in  figure  7.6  and  7.7  show  very  little 

change  in  the  system  gain,  both  are  around  42.0  db,  and  a 
bandwidth  of  about  20  rad/sec.  There  is  a  significant 
change  in  the  phase  diagram  which  is  caused  by  the  zero 
location  shift. 

The  transfer  function  for  Sc  to  depicted  a  gain 

increase  of  about  3  db  for  the  optimized  design  while  the 
phase  remained  similar  for  both  transfer  functions.  In 
transfer  function  Sc.  to  Sa  some  important  aspects  of  the 
problem  are  observed.  The  Bode  diagram  of  the  open-loop 
transfer  function  of  Sc.  to  £&  clearly  indicates  the 
cross-coupling  problems  and  the  pole  placement  and  robust¬ 
ness  design  routine’s  mechanism  of  optimizing  the  system 
gains  to  increase  robustness.  The  gain  is  reduced  from  93 
db  to  about  73  db  and  the  bandwidth  is  reduced  from  above 
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Figure  7.6  Transfer  Function  £&-£q  Design  One,  Nonoptiaized. 

93 


» 


» 


» 


o 


Figure  7.7 


Transfer  Function  Si~Ss  Design  One 


Optiaized. 
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100  rad/sec  to  about  30  rad/sec.  The  reduction  of  bandwidth 
and  gain  in  the  loop  Be.  to  Sb  yields  an  increased  toler¬ 
ance  to  perturbation.  The  transfer  function  for  Sc.  to 
S&  is  shown  in  figures  7.8  and  7.9.  In  the  transfer  func¬ 
tions  for  S&  to  Sc.  the  bandwidth  is  slightly  increased 
from  0.6  to  0.8  rad/sec  and  the  gain  actually  increased  from 
9  db  to  14  db.  The  big  charge  in  the  overall  system, 
however,  is  in  the  transfer  function  from  Sc.  to  Sa  .  This 
is  the  channel  that  the  destabilizing  perturbation  enters 
and  by  greatly  reducing  the  gain  and  bandwidth  in  this 
channel  through  a  change  in  feedback  gain,  the  optimizer 
routine  has  brought  the  entire  system  gains  to  more  balanced 
conditions  and  recovered  a  highly  robust  design. 

The  gain  changes  associated  with  the  robustness  improve¬ 
ment  cause  the  zeros  of  the  various  closed-loop  pole-zero 
diagram  tc  move.  A  comparison  of  the  eight  pcle-zero 
diagrams  is  shown  in  figures  7.10  to  7.13.  The  significant 
feature  of  these  pcle-zero  diagrams  is  the  shift  of  the 
zeros  of  the  optimized  design  in  a  direction  that  attempts 
to  egualize  or  balance  the  frequency  response  for  freguen- 
cies  in  the  vicinity  of  the  minimum  singular  values.  The 
pole-zero  diagram  of  Sb  to  v  will  be  discussed  as  an  example 
of  this  effect.  In  7.10  the  nonoptimized  zeros  are  located 
about  -2  ±  4. 5 j  and  -19.8.  Rhen  the  pole  placement  and 

robustness  routine  has  completed  the  feedback  gain  modifica¬ 
tion  these  zeros  have  shifted  near  -11  and  -6±3j  .  The 
effect  of  these  zero  shifts  is  to  combine  with  the  pole 
locations  tc  egualize  the  frequency  response  as  depicted  in 
figures  7.14  and  7.15.  Zero  shifts  for  the  remainder  of  the 
transfer  functions  provide  similar  results  in  the  ether 
channels.  3y  moving  toward  the  frequencies  associated  with 
the  minimum  singular  values  the  zeros  have  balanced  the 
overall  frequency  response  of  the  system  in  each  channel. 
While  the  channel  gain  modification  is  the  primary  mechanism 
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figure  7,9  Transfer  Punction  Sc  “ Sg  Design  One,  Optimized. 
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figure  7.14  The  $3  to  ▼  Freguency  Hesponse,  Nonopt.. 
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Figure  7.15  The  fig  to  v  Frequency  Response,  Uptiniaed. 
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for  robustness  recovery  the  zero  shift  associated  with  the 
feedback  gain  changes  is  directly  related  to  the  overall 
freguency  response  of  the  system. 

In  design  1  case  2  the  minimum  singular  value  level  was 
increased  from  0.6  to  1.0.  The  pole  placement  and  robust* 
ness  design  procedure  also  gave  a  better  design  for  case  2 

T IBIS  4 

Helicopter  Problem  Feedback  Gains 


by  modifying  the  feedback  gains  as  shown  in  taole  4.  The 
change  in  the  transfer  functions  for  Sc  to  S3  (figure 
7.16  and  7.17)  is  similar  to  the  change  seen  ia  figures7.8 
and  7.9.  The  gain  was  reduced  from  above  90  db  to  about  75 
db  and  the  bandwidth  cut  from  above  100  rad/sec  to  about  25 
rad/sec.  The  added  requirement  of  increased  robustness  did 
not  result  in  a  significant  change  in  performance. 
[Bef.  16]. 
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resign  two  which  is  the  design  that  is  subject  to  param¬ 
eter  uncertainty  in  the  dj?  to  p  channel  was  also  studied 
using  the  pole  placement  and  robustness  design  routine.  In 
this  design  the  singular  values  were  very  low  until  freguen- 
cies  above  100  rad/sec  were  reached.  To  improve  the  design 
it  was  assumed  that  the  pole  locations  that  corresponded  to 
[Hef.  16]  were  the  required  pole  locations  and  that  a 
robustness  singular  value  level  of  0.6  would  provide 
adequate  gain  and  phase  margin  for  the  design.  The  pole 
placement  and  robustness  routine  adjusted  the  gains  in  this 
problem  until  a  minimum  singular  value  of  0.6  was  obtained. 
During  this  adjustment  the  gains  in  channel  $g  were  consid¬ 
erably  reduced  to  offset  the  cross-coupling  between  the  two 
channels.  Again,  in  this  design  as  before  the  channel 
cross-coupling  between  £c.  and  shows  a  marked  change 
in  bandwidth  and  gait  from  above  130  to  35  rad/sec  and  about 
78  db  to  22  db  respectively.  Figures  7.18  and  7.  19  shew  the 
transfer  function  plots  for  this  term.  Figure  7.20  shows 
the  singular  value  improvement.  In  figure  7.21  the  time 
response  is  plotted.  As  can  be  seen  in  the  plot  the 
improvement  in  robustness  for  this  problem  results  in  very 
sluggish  response  and  degraded  performance. 

To  summarize  for  this  problem,  a  given  performance  level 
has  been  chosen  in  terms  of  pole  locations.  The  level  of 
robustness  has  been  set  for  a  desired  gain  and  phase  margin 
based  on  the  universal  gain  and  phase  margin  curve.  The  pole 
placement  and  robustness  routine  has  been  able  to  improve 
the  robustness  level  by  changing  the  feedback  gains  that 
affect  the  channel  cross-coupling.  This  robustness 
recovery  is  affected  by  modification  of  the  system  feedback 
gains  in  such  a  manner  that  cross  coupling  gains  are  reduced 
so  that  small  cross-coupling  perturbations  do  not  drive  the 
system  into  instability.  The  open-loop  transfer  function 
plots  have  been  used  to  indicate  how  this  mechanism  operates 
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Pigure  7.18  Transfer  Function  $c-  $8  Design  Two,  Nonopt.. 
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Figure  7-20  Singular  Value  Plot  Design  Two. 
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and  have  teen  shown  tc  be  an  alternative  indicator  of  chan¬ 
nels  that  may  be  affected  by  cross-feed  perturbations.  The 
pole-zero  diagrams  of  the  closed-loop  transfer  functions  of 
the  transfer  matrix  further  indicate  that  zero  movement  is 
in  a  direction  that  equalizes  the  gain  level  of  the 
frequency  response  curves  in  the  vicinity  of  the  lowest 
singular  values  providing  a  more  balanced  system  response. 


fill.  SI|£L£  OBSERVER 


The  pole  placement  and  robustness  design  procedure  can 
also  be  used  for  robustness  recovery  in  observer  design. 


Figure  6.1  Simple  Observer. 


Given  an  observer  as  represented  in  figure  8.1  it  has  been 
shown  that  the  system  differential  equation  may  be  written 
as  equation  8.1 

i  «fi  -BF]  jV)  ♦  [b]  r  (8.1) 

$  [kC  Ac-BFj  Jij  [3 

At  =  A  -  KC 

where  x  is  the  state  and  x  is  the  observer  variaole,  and  the 
solution  may  be  separated  to  independent  solutions  for  the 
feedback  gains  and  the  observer  gains.  This  separation 
allows  the  use  of  the  feedback  gains  to  set  the  pole 
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locations  and  then  a  second  optimization  run  with  set  feed- 
tack  gains  so  that  the  observer  gains  may  be  computed  to 
adjust  tie  system  rotustness  level. The  observer  pole  loca¬ 
tion  could  also  be  placed  using  the  observer  gains,  K.  In 
the  currently  implementation  of  the  pole  placement  and 
robustness  design  routine  the  observer  poles  are  simply 
restricted  to  areas  cf  the  left  half  plane. 

In  this  chapter  a  simple  stable  observer  system  will  be 
analyzed  based  on  [Bef.  17].  Given  the  system,  equation 
8.2, 


5  -f°  ’l  5  +  [°1  u  *  f  351 

t3  -«J  LlJ  1-6  ij 


(8.2) 


y*  [2  1]  x  + 

where  E  fy)  *E  =0  and  E  (jU  )j  (f) )  -E  (t )  f  (r) )  =  £ ( t-r) 
the  feedback  law  of  equation  8.3 


with 


u  *  -  (5C  10)  x  ♦  50r 


(8.3) 


An  analysis  has  been  done-  to  compare  results  of  the  numer¬ 
ical  optimization  procedure  with  results  presented  in 
[Bef.  17].  The  results  for  an  optimal  regulator  design 
using  quadratic  cost  criteria  as  stated  in  equation  8.4  are 
given  in  table  5 


=  /uV 


Hx+  u8)  dt 


(8.4) 


with  h  =  4  y?  (  y55  1) 

Figure  8.2  shows  a  Nyquist  plot  of  the  full  stats  regulator, 
the  optimal  filter  and  a  fast  filter.  The  full  state  design 
had  poles  at  s=  -7.0±j  2.0  and  feedback  gains  of  50  and  10. 
The  optimal  filter  as  shown  in  table  5  had  poles  of  -7.0  ±  j 
2.0  with  gain  and  phase  margins  of  -6.75  db  and  115  degrees. 
The  optimum  filter  gains  were  30  and  -50.  Using  a  faster 


Figure  8.2  Byguist  Plot. 

filter  also  gives  poor  gain  and  phase  Margins.  The  gain 
nargins  are  on  the  order  of  -.98  while  the  phase  aargin  is 
less  than  10  degrees.  The  bandwidth  also  increased  from  12 
to  40  rad/sec.  k  recovery  procedure  based  on  a  aodif ication 
to  the  process  noise  aatrix  [Bef.  17]  may  be  applied  to  the 
problem.  This  procedure  can  recover  a  large  amount  of  the 
robustness  that  was  lest  with  the  observer  addition.  Figure 
8.3  shows  data  obtained  for  several  trials  of  the  fictitious 
noise  procedure.  The  gain,  phase  margin  and  ether 
parameters  may  be  found  in  table  5. 


TABLE  5 

Observer  Parameter  Data 
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The  singular  values  o£  the  state  feedback  and  optiaal 
ebserver  systems  were  computed  for  comparison  # ith  results 
produced  using  the  pole  placement  and  robustness  design 
recovery  procedure.  Figure  8.4  indicates  a  loss  in  robust¬ 
ness  represented  by  the  observer  singular  values.  lover 
singular  values  are  less  robust.  For  this  single- input 
single-output  observer  the  Hyguist  diagram  will  be  used  to 
define  gain  and  phase  margins.  Since  a  singular  value  of  1.0 
is  indicative  of  a  linear  guadratic  level  of  gain  and  phase 
margin, i.e.  GM=-6  db,  and  PN=±60  degrees,  this  was  chosen 


figure  8.3  lyguist  Plot  for  Robustness  Recovery. 

as  the  design  level.  The  pole  placenent  and  robustness 
routine  was  used  to  recover  robustness  while  setting  the 
pole  placenent  at  the  state  feedback  pole  locations  of 
-7±j2.  The  nininun  input  singular  value  level  was  set  at 
1.0.  The  pole  placenent  and  robustness  routine  was  also  set 
to  place  the  observer  poles  anywhere  between  -100  and  -2 
that  would  provide  robust  design.  The  plant  poles  were 
placed  at  -7.05  ±j  1.82.  The  feedback  gains  for  this  run 
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SINGULAR  VALUES 


Figure  8.4  Singular  Values  of  Observer  Systeu. 
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were  50  and  10.09.  The  optimizer  output  produced  a  filter 
with  a  gains  of  1.00018  and  -2.02742.  The  observer  pole 
locations  were  -1.99  ±  j  .001,  very  near  the  plant  zero 
location.  A  significant  change  was  produced  in  singular 
values.  Figure  8.5  shows  all  singular  value  curves  plotted 
together.  The  optimizer  solution  for  this  problem  is  well 
above  the  optimal  filter  curve  at  low  frequency.  Figure  8.6 
shows  the  Nyquist  plot  of  the  optimizer  developed  design. 
The  system  has  a  gain  margin  of  -6  db  and  a  phase  margin  of 
117  degrees.  The  most  significant  differences  between  the 
two  designs  being  that  the  observer  poles  are  close  to  the 
plant  zero  locations  and  the  filter  gains  are  much  lower  for 
the  optimizer  solution. 

Using  the  OPTSYS  program  with  the  pole  placement  and 
robustness  routine  computed  gains  as  the  design  parameters 
the  data  for  the  observer  filter  was  computed.  The  error 
covariance  matrix  was  found  to  be 
306.5  -456.5 

-456.5  7C0.1 

These  values  compare  favorably  with  the  trends  established 
in  table  5.  The  last  comparison  of  the  pole  placement  and 
robustness  design  recovery  procedure  was  an  analysis  of  the 
time  response  curves.  Figure  8.7  shows  the  comparison  plot. 
The  design  obtained  using  the  pole  placement  and  robustness 
routine  did  not  degrade  system  performance. 

Cne  additional  analysis  was  conducted  that  set  the 
desired  design  parameters  slightly  differently.  In  this  run 
the  pole  placement  and  robustness  routine  was  set  to  place 
the  observer  poles  between  -10  and  -100.  The  pole  placement 
and  robustness  routine  was  unable  to  totally  satisfy  this 
requirement.  It  violated  one  of  the  constraints  and  moved 
the  larger  observer  pole  to  -2.25  which  is  near  the  optimum 
pole  location  and  also  the  high  g2  values  of  the  fictitious 
noise  procedure.  The  smaller  pole  was  moved  out  to  -70.8 
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SINGULAR  VALUES 


Figure  8.5  Singular  Value  Comparison  Plot. 
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Figure  8.6  Hyguist  for  Coiputed  Robustness  Recovery 


which  corresponds  to  the  same  level  of  pole  movement  found 
when  g*  fceccmes  large.  The  singular  value  levels  were  raised 
significantly  as  shown  in  figure  8.5.  The  state  and  error 
covariance  matrices  were  of  the  same  order  as  the  matrices 
found  in  table  5.  This  design  has  a  phase  margin  cf  58 
degrees.  The  gain  margin  was  about  -19  db.  Even  though  seme 
of  the  constraint  conditions  on  the  design  were  not  met  the 
design  demonstrates  excellent  robustness  recover/  . 


AMPLITUDE 


The  pcle  placement  and  robustness  design  routine  has 
produced  a  robust  design  for  this  observer  based  system. 
This  design  is  obtained  by  using  a  numerical  optimization  to 
directly  manipulate  the  feedback  and  filter  gains. 
Modification  of  the  IQ  functional  eguation  as  done  in  the 
fictitious  noise  procedure  is  not  required.  The  pole  place¬ 
ment  and  robustness  routine  solution  for  this  problem  has 
lower  gains  than  those  found  using  the  fictitious  noise 
adjustment.  The  routine  provides  a  good,  direct  methodology 
for  selecting  the  feedback  and  filter  gains  for  a  rcbust 
observer  design  with  excellent  performance. 


a 

m 
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IX.  BOBOSI  OBSESTER  DESIGN 

This  chapter  will  be  devoted  to  a  short  discussion  of 
the  robustness  recovery  of  a  fourth  order  observer  based 
controller.  The  problem  is  a  helicopter  problem,  £Bef.  6], 
In  this  case  the  helicopter  model  is  that  of  the  longitu¬ 
dinal  control  loop  of  a  CH-47.  The  nominal  model  is  taken  to 
be  the  system  of  equations  9.1  and  9.2  for  an  aircraft  speed 
of  forty  knots. 


• 

X 
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-32. 

x  + 

0.  14 
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-0.  14 
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-30. 

.36 

-8.6 

0. 

0.018 

-1.6 

1.2 

.35 

.009 

0. 

m 

0. 

1. 

0. 

J 

0. 

0. 

y  =  C  1  0  0  ]  x  (9.2) 

0  0  0  57.3^ 

In  this  problem  the  controller  is  formulated  as  in  figure 
9.1  which  leads  to  an  open-loop  transfer  function  of  the 
form  of  eguation  9.3 

K  (S)G(s)  =F  (sI-A+  EI+KC)  "  1  RC  (si- A)  " 1  B  (9.3) 

The  pole  placement  and  robustness  recovery  procedure  was 
applied  to  this  problem. 

First,  the  standard  full  state  feedback  design  was 
carried  out  using  the  Naval  Postgraduate  version  of  OPTSYS. 
This  design  produced  excellent  singular  value  output  for  the 
return  difference  as  shown  in  figure  9.2.  The  lowest 
singular  value  being  essentially  1,  corresponding  to  a  LQ 
design  with  -6  db  to  CD  gain  margin  and  60  degrees  of  phase 
margin.  The  time  response  of  the  system  was  good  as  shown  in 
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Figure  9.1  Observer  Based  Controller. 

figure  9.3,  reaching  steady-state  in. about  four  seconds  with 
only  a  slight  overshoot.  Since  the  characteristics  for  this 
system  are  acceptable  no  further  design  iteration  was 
carried  out.  The  full  state  feedback  became  the  baseline 
design.  Assuming  that  full  state  feedback  vas  not  available 
and  only  two  measurements  could  be  produced  an  observer  vas 
developed  to  control  the  two  measured  outputs,  vertical 
velocity  and  pitch  attitude.  Using  the  measurement  matrix  of 
eguation  9.2,  OPTSYS  vas  used  to  develop  an  optimal  observer 
system  for  this  problem.  The  singular  values  of  the  return 
difference  matrix  and  the  time  response  were  plotted  for 
comparison  to  the  full  state  design.  These  are  shovn  in 
figures  9.2  and  9.3  respectively.  Hote  that  the  singular 
value  of  the  return  difference  matrix  is  as  lov  as  0.16  at  U 
rad/sec.  This  equates  to  a  gain  margin  of  -1  db  to  1  dt  and 
less  than  10  degrees  phase  margin.  The  time  response  is 


P 
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NOTE:  ALL  RESPONSES 

:  ESSENTIALLY!  IDENTICAL 
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TIME  -  SEC 


Flotted  in  figure  9.3.  The  steady-state  was  reached  for  the 
velocity  after  a  slight  overshoot  in  about  four  seconds. 

As  the  final  step  in  the  analysis  the  pole  placement  and 
robustness  routine  was  employed  to  recover  robustness  of  the 
observer  based  system.  The  pole  placement  and  robustness 
routine  was  first  used  to  set  the  poles  at  approximately  the 
same  location  as  the  poles  of  the  LQ  regulator.  These  were 
assumed  to  be  the  desired  pole  locations.  The  pole  placement 
and  robustness  routine  was  then  used  to  vary  the  filter 
gains,  K,  until  the  desired  level  of  robustness  was  reached. 
The  desired  singular  value  level  chosen  was  1.0  which  corre¬ 
sponds  tc  -6  db  to  eo  db  gain  margin  and  60  degrees  of  phase 
margin.  Figure  9.2  shows  that  the  pole  placement  and  robust¬ 
ness  procedure  failed  to  totally  recover  the  robustness 
level  to  1.0.  The  minimum  singular  value  reacaed  was  only 
C.936.  This  singular  value  equates  to  a  gain  margin  of  -5.0 
db  to  eo  and  a  phase  margin  of  55  degrees.  While  this  is 
slightly  less  than  the  design  objective,  it  is  far  superior 
to  the  optimal  observer  design  discussed  in  the  previous 
paragraphs.  This  recovery  was  made  by  making  an  optimiza¬ 
tion  run,  finding  the  filter  gains  with  near  zero  values  and 
freezing  these  values  at  zero.  This  reduced  the  number  of 
design  variables  the  pole  placement  and  robustness  routine 
was  required  to  manipulate  in  a  second  optimization  run  and 
gave  a  higher  robustness  solution. 

This  analysis  and  the  second  order  observer  analysis 
presented  earlier  clearly  indicate  that  it  is  possible  to 
use  the  pcle  placement  and  robustness  procedure  and  separa¬ 
tion  principle  to  develop  a  robustness  recovery  procedure. 
The  pcle  placement  and  robustness  routine  provides  robust¬ 
ness  recovery  by  direct  modification  of  the  feedback  and 
filter  gains.  This  procedure  requires  no  modification  to  LQ 
cost  functionals  or  ether  parameters  as  done  in  the  ficti¬ 
tious  noise  adjustment  method  commonly  used  for  robustness 
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recovery.  By  providing  direct  gain  adjustment  the  pole 
placement  and  robustness  procedure  results  in  a  practical 
design  with  relatively  low  observer  gains  and  good  perform¬ 
ance.  The  procedure  is  simple  and  straight  forward  with  the 
only  difficulty  being  a  requirement  to  sometimes  modify 
initial  starting  values  or  optimizer  codes  to  force  the 
solution  toward  the  desired  point. 


X.  C0NC1PSI0NS 


An  effective  method  of  robustness  multivariable  control 
design  utilizing  a  numerical  optimization  based  algorithm 
has  teen  developed.  The  pole  placement  and  robustness  design 
routine  coupled  with  the  Automated  Design  Synthesis  program 
provides  the  designer  an  excellent  tool  with  which  to  attack 
the  robust  design  problem. 

The  pole  placement  and  robustness  design  routine  has 
demonstrated  the  capability  of  providing  designs  that  solve 
the  problems  caused  by  cross-coupling  perturbations  which 
reduce  robustness  in  multivariable  systems.  This  design 
improvement  is  accomplished  by  modifying  the  system  feedback 
gains  in  such  a  manner  that  the  gain  in  channels  that  are 
affected  by  cross-coupling  perturbations  is  equalized  with 
other  system  gains  to  reduce  this  cross-coupling  effect.  The 
gain  changes  are  accompanied  by  zero  shifts  which  also 
influence  the  gain  distribution  and  freguency  -response  of 
the  system. 

Perturbation  problems  in  multivariable  systems  have  been 
shown  to  be  detectable  by  singular  value  analysis  and  by 
using  the  Bode  magnitude  diagram  of  the  open- loop  transfer 
functions  of  the  system.  In  the  open  loop  transfer  function 
large  differentials  in  Bode  gains  and  bandwidths  are  indica¬ 
tive  of  problem  areas  for  cross-coupling  perturbations. 
Eobustness  is  obtained  by  the  pole  placement  and  robustness 
design  program  by  modifying  those  gains  and  bandwidths  asso¬ 
ciated  with  the  cross-coupling  perturbations  thus  reducing 
the  amount  of  energy  coupled  from  the  perturbation  into 
ether  channels.  An  associated  zero  shift  has  been  observed 
when  these  gain  modifications  take  place.  This  zero  shift  is 
in  the  direction  of  poles  that  are  located  in  the  vicinity 


130 


cf  the  frequency  of  the  minimum  singular  value  and  tends  to 
equalize  the  frequency  response  curve  gains  in  this  region. 

The  use  of  numerical  optimization  to  recover  robustness 
in  observer  based  designs  was  demonstrated.  The  pole  place¬ 
ment  and  robustness  routine  was  applied  to  problems 
previously  solved  using  the  fictitious  noise  procedure  for 
robustness  recovery.  The  direct  manipulation  of  feedback  and 
filter  gains  by  the  pole  placement  and  robustness  routine 
provided  a  highly  robust  design  with  relatively  lew  filter 
gains.  The  problem  of  robustness  recovery  in  filter- 
observer  designs  has  been  solved  in  a  straight  forward  and 
highly  practical  manner. 
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